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Abstract 

A Sensorless Backstepping technique combined with a 
field orientation scheme has been developed for the 
control of induction motor to achieve rotor angular speed 
and rotor flux amplitude tracking objectives. An accurate 
knowledge of the rotor speed is the key factor in 
obtaining a high-performance and high-efficiency 
induction-motor drive. For that, a Model reference 
adaptive system (MRAS) rotor-speed estimator has been 
incorporated for which stability and robustness are 
assured. Simulation results are presented to validate the 
effectiveness of the overall proposed sensorless control 
scheme under the extreme conditions. 



Keywords: Induction motor, Backstepping design, Vector 
control, MRAS speed identification. 



Nomenclature 



r x ,r 2 


stator, rotor resis tan ce 


h , h 


stator, rotor current 


X x ,X 2 


stator, rotor flux linkage 


^2 d 


amplitude of rotor flux linkage 


u x ,u 2 


stator, rotor voltage input 


CD 


rotor angular speed 


p 


angle between the rotor flux linkages 


n p 


number of pole pairs 


L x ,L 2 , L m stator , rotor and mutual inductance 


T v T 2 


stator, rotor time cons tan t 



J ,T l inertia, load torque 
(A0 9 (•) in (d,q) frame 
(■Ui (•)'« {a, b) frame 
(?) estimate of (.) 



o = 1 - - 



h 2 = 



L x L 2 



2 

gL, L 2 



a = - 



p= 






^ n p L m 

2JL-, 



gL, L 2 
d/ 



hi * 



A. 

gLi 



>P=** 



/ dt 



1. Introduction 

Induction motors have been widely used in high 
performance alternating current (ac) machine drives, 
requiring speed information. Introducing shaft speed 
sensor decreases system reliability. Different solutions 
for speed sensorless control of the induction motor (IM) 
have been proposed [1, 2, 3]. Among them, MRAS based 
techniques have been proven to be one of the best 
methods being proposed by the researchers, due to its 
design and implementation simplicity and high 
performance capability [4]. 

In the last two decades, many modified nonlinear state 
feedback schemes such as input-output feedback 
linearization [5], passivity-based control [6, 7] and 
Sliding-Mode (SM) control [8] have been applied to the 
IM drive. 

Specially in the last few years, in the field of adaptive 
and robust control, there has been a tremendous amount 
of activity on a special control scheme known as 
“Backstepping” [9, 10, 11]. With adaptive Backstepping 
design technique researchers were able to successfully 
design controllers which achieved global stabilization in 
the presence of uncertain parameters[12, 13]. Using this 
design methodology the construction of both the 
feedback control laws and associated Lyapunov functions 
is systematic. Strong properties of global or local stability 
are built into the nonlinear system with uncertain 
parameters [14]. 

In this paper, we have integrated a powerful tool of 
control (Backstepping) with a powerful tool of speed 
estimation (MRAS). To validate the effectiveness of this 
integration, the overall scheme is subjected to the 
extreme variation of speed reference and the benchmark 
test. 

The rest of this paper is organized as follows: in section- 

2, we present the IM model used and the Backstepping 
control design, in section-3, we present the rotor flux 
estimator, in section-4 the MRAS speed estimator is 
detailed and in section- 5 the stability of the overall 
scheme is discussed. 
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2. Model Description and Control Design 

Assuming linear magnetic circuits and balanced three 
phase windings, the fifth-order nonlinear model of IM, 
expressed in the stator frame is [15]: 



dco 3 n p L m , . . \ T l 

“77-^77 v^lahb ~Hbha) T 

at 2JL2 J 



dX 



dt 

dX 



A 

Li 



— = ~—X 2 a - n p a>X 2b + — L„, i 



L 2 



'm L \ a 



-TT = ~^2b + n p^2a + Y^ L mhb 
at L 2 y L 2 



dha _ k m^2 , , n p^m , 

7T k 2a + _ r 7 {i>k 2 b ' 



dt gL x L 2 



oL^L 2 



1 IR 2+&1 

al^Ll 



U a 



crfy 



Ha 



di\b _ L m R 2 7 n pL m * L m R 2 + L 2 R x 

~7 ~ ~ ~ t 2 ~^2a 1 

at aL, L 2 gL| L 2 aL, L 2 



lb 



1 

crfy 



(i) 



Hb 



We can see that (1) is highly coupled, multivariable 
and nonlinear system. It is very difficult to control the IM 
directly based on this model. State transformation to 
simplify the system representation is required. A well- 
known method to this end is the transformation of the 
field orientation principle. It involves basically a change 
of the representations of the state vector {i la d xb ,X 2a ,X 2b )m 
the fixed stator frame (a,b) into a new state vector in a 
frame (d,q) which rotates along with the flux vector 
(Aa’Ab) • 



Mathematically, the field oriented transformation can 
be described as: 



• = ^2 aha + ^2bkb • 

4 2 d = d%a + %b> H, = P = arctan 



^lahb ~ ™ Ibha 



Va 



+ X- 



'2 b 



^2 b 

^2 fly 



( 2 ) 



Where much simplification is gained by the fact 
that X 2q = 0 . 

Using this transformation and the notations in the 
Nomenclature, the state equations (1) can be rewritten in 
the new state variables as: 



<Zco 

dt 



~ M^ 2 dhq " L 



Tl 

J 



di i 



-^- = - r lih ? ~$n p (»X 2d -n p a>i ld -R 2 



r\ 2 I u .+aL n 



l \q l \d 
^2 d y 






dk 



ld = -aR 2 X 2d + aL m R 2 i u 



dt 

di ld . . „ 

~ _r li hd + n p™ l \q + R 2 



-n2 I ld+ a ^2d+ aL mX L - 



1 

afy 



Hd 



— = no + aL R 2 — 

dt p m 2 x 2d 



2d 



(3) 



From this structure, the Backstepping control 
mechanisms for the rotor angular speed regulation and 
the flux generation can be better applied to replace the 
traditional nonlinear feedback plus (proportional-integral) 
PI control of the field oriented control technique for 
better performance. 

The basic idea of the Backstepping design is the use of 
the so-called “virtual control” to systematically 

decompose a complex nonlinear control design problem 
into simpler, smaller ones. Roughly speaking, 

Backstepping design is divided into various design steps. 
In each step we essentially deal with an easier, single- 
input- single-output design problem, and each step 
provides a reference for the next design step. The overall 
stability and performance are achieved by a Lyapunov 
function for the whole system [16]. 

Now we turn to the Backstepping design steps. First, for 
the rotor angular speed, and rotor flux amplitude tracking 
objectives, define the tracking errors as 

e \ = 03 ref ~ 03 ( 4 ) 

^3 — ^ ref ^2 d 

Then the error dynamical equations are 

T 

k =(°ref-^2dh q +-7- 

J (5) 

e 3 = X ref + aR 2 X 2d -a L m R 2 I u 



Since our objectives require that the two errors converge 
to zero, we could satisfy the objectives by viewing X 2d i lq 

and R 2 i ld as “virtual control” in the above equations and 
use them to control e x ,e 2 . Certainly, X 2d i ]q is only a 

constant from the torque generated by the motor, and we 
include R 2 in the virtual control R 2 i ld . We use the 
Lyapunov function 



V = -el + -el 

2 1 2 3 



( 6 ) 
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To derive stabilizing feedback virtual controls A 2 d i ]q and 
R 2 i ld as follows. The derivatives of V along the error 
equations ( 5 ) is 



V = e x e x + e 3 e 3 



= e, 



co 



ref 



T 

V-^ldhq + _ T 



+ e 3 [i ref + aR 2 X 2d - a L m R 2 i ld ] 



= -k x e x -k 3 e 3 + e x 



+ e 



\l c 3 e 3 + 



T 

k x e x +a r ef-\^2dhq+-j- 



( 7 ) 



ref + aR 2^2 d ~ aR m R 2hd 1 



Where k { ,k 3 are positive design constants that determine 

the closed-loop dynamics. Thus the tracking objectives 
will be satisfied if we choose 



( 8 ) 



/ . \ 1 


T,~ 


\ 2 dhq) re j2 ^ 


K e \ + & ref + — 



( 9 ) 



( R 2 hd )ref ~ ~~T~ i k 2 e 3 + ^ ref + a R{k ld ] 
CCA,,, 



since the choice of (8) will gives 

V = ~k x e\ - k 3 e\ < 0 

The equations (8) indicate what the virtual controls 
should be that in order to satisfy our control objectives. 
So they provide references for the next step of the 
Backstepping design, which is essentially to make the 
"torque” and i u behave in the desired manner. 

Now we go to next step and try to make the signals A 2 d i ]q 
and i u behave as desired. So we define again error 
signals involving the desired variables in (8) 

e 2 - {k 2d h q ) ref ~^2 dhq 



~ ~ [k\ e l + ®ref ] + “T “ dh 

p \U 

“ ( R 2 Hd )ref ~ R 2 hd 



4 



a L 



7 — [ k 3 e 3 + k ref + aR 2 k 2d ]~ R 2 hd 



( 10 ) 



Then the error equations ( 5 ) can be expressed as 
<?i = -Vi + H e 2 

e 3 = -k 3 £ 3 + aL m e 4 (n) 

Also, the dynamical equations for the error signals 



2 ,e 4 can be computed as 



1 



^2 “ ^1 ~~^j^2 d U \q 

^4 = ^2 — R 2 U \ d 

GL^ 



( 12 ) 



Where ' s(i = 1 , 2 ) are known signals that can be used 
in the control, and their expressions are as follows 



<h =— (- Vi +ve 2 )+-<»ref + r llWl 9 

0 n 

+ fin p a\ 2 2d + n p aX 2 d i ld +/?2( a + r l2)'i ? ^2rf ( 13 ) 



<l>2 =-T _ (- fe 3 e 3 + aL m e 4)+-^— k 



aL„. 

- R 2 {-nihd+ n p a hq)- R : 



a 

f 



ref ~ (^ 2 d ~ L mhd ) 

An 

•2 A 



-r\ 2 i ul +afA 2 d +aL in -^- 

k 2d 



Now we extend the Lyapunov function in (6) to include 
the state-variables e 2 , e 4 as 



K - — [ e i + e 2 + e\ + e\ 



( 14 ) 



We use V e to derive both the control algorithm. To this 
end, we again compute the derivate of V e along with the 
error equations (1 1), (12) 

K = e A + *2*2 + *3*3 + *4*4 

= e\(-k l e ] + \ie 2 ) 



+ 



. ^2 d U 



gL^ 






+ e 3 (-k 2 e 3 +aL m e 4 ) 

/ ' 

1 



+ £4 



^2 



gI^ 



R 2 u \d 



( 15 ) 



— k x e x ^2^2 k^3 ^4^4 



+ ^3 



+ 



(a e x +k 2 e 2 + f —X 2 d u 1 

GZ4 



aL m e 2 +k A e A +§ 2 —R 2 u ld 

gl^ 



From the above we can obtain the control laws as 

= ^ L t le l+M2+ ( l>l] 



4 lq ' 



' J 2d 



u \d = ^f[ aL m e 3 + k 4 e 4 + <t> 2 ] 



R 



Witch leads to 



V e — k\&\ & 2*2 k 3 (? 3 k 4 e 4 ^ 0 



( 16 ) 



( 17 ) 



3. Adaptive Rotor Flux Estimator 

The rotor flux has been estimated using the simplified IM 
model ( 3 ) obtained after the application of the field 
orientation principle transformation. Where much 
simplification is gained by the fact that X 2q = 0 . This 

allows estimation of the rotor-flux vector in polar form as 
given by 



^ 2 d ~ ~ aR 2 ^ 2 d + aR m R 2 hd 



( 18 ) 



4. MRAS Speed Identification 

It is based on the MRAS identification approach [ 17 ]. 
Since the motor voltages and currents are measured in a 
stationary frame of reference, it is also convenient to 
express these equations in the stationary frame (a,b) [2]: 
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^2 a 



v 2 b 



A 



*1 a 



( R { + OLy p ) 0 

0 ( R x + gZ^ p ) 



l \a 



b_ 



^2 a 




1 

3 

1 

E? 

\ 

1— H 
1 

1 


^2 a 




Ha 


^2 b _ 




03 (~l/T 2 ) 


<N 


T 2 


}\b _ 



(19) 

(20) 



Figure 1 illustrates the estimation of the motor speed by 
means of MRAS techniques. Two independent observers 
are constructed to estimate the components of the rotor 
flux vector: one based on (19) and the other based on 
(20). Since (19) does not involve the rotor speed to be 
estimated, this observer may be considered as a reference 
model. 

Within this approach an identifier consists of a tuning 
model (20) depending on an estimate of the unknown 
rotor speed and a mechanism to adjust the estimate. This 
adjustment is performed to make the output of the tuning 
model asymptotically match the output of the reference 
model. 




Figure 1: Structure of MRAS speed estimator 



K p and K t : are the gain constants of the adaptation 
mechanism. 

5. Stability Analysis 

To prove the convergence properties and stability of the 
overall closed loop system, we will prove that all the 
signals of the states are bounded and all tracking errors 
and all estimation and observation errors converge to 
zero asymptotically. For that, the proof is established for 
both adaptive Backstepping controller and adaptive 
MRAS identifier. 

a. Adaptive controller: 

The proof is established via Barbalaf s Lemma using the 
Lyapunov function candidate [16]. 

For that, we have to prove 

e 4 £ £oo n £a, h , ^2, £ 3 * £4 € £00 (22) 

These conditions imply, via Barbalaf s Lemma [10] that 
the closed-loop is stable and that e x ,e 2 ,e 3 and e 4 all 
converge to zero. 

We begin with the facts that V e is positive definite, and 
V e is semi-negative definite, thus all the error signals 
e\, e 2y e 3 and e 4 are pounded. With our assumptions on the 
reference signals co ref , ,A ref ,we get that the signals 4>i and 
4>2 are bounded. Therefore by equations (11), (12), we 
know that e x ,e 2 ,e 3 and e 4 are bounded. Now, what left is 
to prove 

«li«2 1*8164 6 £2 



The adaptive mechanism should be designed to ensure 
the stability of the control system. 

Rotor speed is obtained from the adaptation mechanism 
as follows: 



00 = 



K p + 



K, 



(^2 /A 



2 a ^2a^2b j 



( 21 ) 



Figure2 shows the block diagram of the speed 
identification MRAS based on this adaptation mechanism. 




Figure2: Block diagram of MRAS speed estimator 



By the equation (17) we get 

V e ~ - kie\ - k 2 e 2 - < -k\e\ (24) 

Integrating on both sides from 0 to 00 we have 



V e (oo) - V' c (0) < - / hie\ dt 
Jo 

US 

V*(G) -Ve(oc) >_ he* 
Jo 



dt 



(25) 



(26) 



which imply that e x In the same way, we can also 
prove that e 2 , e 3 , e 4 &£ 2 - 

Therefore, we have concluded that all the signals of the 
states are bounded and that the error signals e x , e 2 , , e 3 , e A 
converge to zero asymptotically. 



b. Adaptive identifier 

For the MRAS identifier design, it is important to take 
account of the overall stability of the system and to 
ensure that the estimated quantity will converge to the 
desired value with suitable dynamic characteristics. 
Landau [18] has described practical synthesis techniques 
for MRAS structures based on the concept of 
hyperstability. When designed according to these rules, 
the state error equations of the MRAS are guaranteed to 
be globally asymptotically stable. 

Subtracting (20) from (19), we obtain the following state 
error equations: 
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d 


e a 




7-1 /T 2 ) -CO 


e a 


_|_ 


i 

<N 

1 


dt 


Jb_ 




to f-i/r 2 ) 


Jb. 


1 


_^2 a 



That is, — [e] = [a][<?] - \\V ] (28) 

dt 




Figure3: Equivalent non-linear feedback system 



The estimation algorithm (adaptation mechanism) is 
chosen according to Popov’s hyperstability theory, 
whereby the transfer function matrix of the linear time 
invariant subsystem must be strictly positive real and the 
nonlinear time varying feedback subsystem satisfies 
Popov’s integral inequality, according to wich 
h 

^(e T W^dt > 0 in the time interval [0,q] for all 
o 



h - 0 [19]. 

Since ob , is a function of the state error, these equations 
describe a nonlinear-feedback system as illustrated in 
figure3. Following Landau, hyperstability is assured, 
provided that the linear time-invariant forward-path 
transfer matrix is strictly positive real and that the 
nonlinear feedback (which includes the adaptation 
mechanism) satisfies Popov’s criterion for hyperstability. 
Popov’s criterion requires a finite negative limit on the 
input/output inner product of the nonlinear feedback 
system. Satisfying this criterion leads to a candidate 
adaptation mechanism as follows: 

Let cb = 0 2 ([^])+ f ®i([e])iT (29) 

Jo 

Popov’s criterion requires that 

f > -/q for all t x > 0 (30) 

Jo 



where yl is a positive constant. Substituting for \e\ and 
\\V ] in this inequality and using the definition of cd , 
Popov’s criterion for the present system becomes 



r 

JO 



\ e a^2b ' 



' e b^2a 



] ®-0 2 ([e])-J 

JO 



■dt > -y 0 
(31) 



Using this expression, it can be shown that Popov’s 
inequality is satisfied by the following functions: 



(e b A 2a -e a ^2b)~ {^2b^2a ~ ^2 a^2b) 

®2 “ K p (e b A 2a ~ e a^2b)~ K p {^2b^2a ~ ^2a^2b ) 



(32) 



6. Simulation Results 

The overall configuration of the control system for IM is 
shown in Figure4. The effectiveness of the proposed 
controller combined with the rotor speed estimation has 
been verified by simulations in Matlab/Simulink. 

The parameters of the induction motor used are given in 
Appendix. The simulation results have been obtained 
under a constant load of lON.m between (2 and 6)s. 
Parameters of the Backstepping controller are: 
kl= 28000, k2= 12000 , k3= 21000 and k4= 13000. 
Parameters of the MRAS speed identifier 
are: K p = 50000 , K t = 100000 . 

The simulation results of the performance tracking are 
demonstrated in the following cases: 

Case 1. The speed references are smooth functions, 

G>ref = 50(l -exp(-2;tf))+10sin(2;zt) and 50sin(2;tf) 

The reference flux is set to 0.5 wb. 

As shown in the results of Figure5, the performance are 
satisfactory. 

Case 2. The speed reference is a sort of step type (200 
rad/sec) to validate the proposed controller. 

The reference flux is set to 0.5 wb. 

Figure 6 shows the performance of the proposed scheme 
with a step speed reference. The errors of speed 
estimation and speed tracking are rapidly converged. 

The control objective is achieved even with the crutial 
condition. 

Case 3. The speed tracking controller is operated in a 
critical situation of benchmark commands (rapidly 
changes as 30 - 0 — 30 - 0 - 10 rad/sec). 

Figure 7 shows the satisfactory performances of the 
speed tracking. We can see that the actual speed follows 
the speed command. Although the speed reference 
rapidly changes, the actual speed output of the induction 
motor is closely follows up. On the other hand, the 
current input is zero as the motor stops. 

Thus, the simulation results confirm the robustness of the 
proposed scheme with respect to the variation of the 
speed rotor reference. 

7. Conclusion 

In this paper, a novel scheme for speed and flux control 
of induction motor using online estimation of the rotor 
speed has been described. The nonlinear controller 
presented provides voltage inputs on the basis of only 
stator currents measurements and guarantees rapid 
tracking of smooth speed and rotor flux references. 
MRAS rotor speed estimator has been used. In simulation 
results, we have shown that the proposed nonlinear 
adaptive control algorithm achieved desirable tracking 
performance within a wide range of the operation of the 
IM. The proposed method also presented robustness 
properties with respect to the extreme variation of the 
rotor speed reference. The other interesting features of 
the proposed method are that it is simple and is easy to 
implement in real time. 

It will be noted that the presence of the pure integrators 
in MRAS scheme brings the problems of initial 
conditions and drift. A low pass filter was used to replace 
the pure integrator [2]. 
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APPENDIX 

INDUCTION MOTOR DATA 



Stator resistance 
Rotor resistance 
Mutual inductance 
Rotor inductance 
Stator inductance 
Number of pole pairs 
Motor load inertia 



1.34 ohms; 

1 .24 ohms; 
0.17 H; 

0.18 H; 

0.18 H; 

2 

0.0153 kgm 2 ; 
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Figure4: Overall bloc diagram of the Sensorless control scheme for IM 





Time (sec) Time (sec) 





Figure5: Tracking performance 



Figure6: Tracking performance Figure7: Tracking performance 



Left : o) re f = 50(l — exp(- 2m)) +10 sin(2^r) 
Right : a) rpf =50sin(2;tf) 



Step type speed tracking 
co re f = 200 rad /sec 



Benchmark speed tracking 
v re f =(30 0 -30 0 10)ra<i/sec 
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Abstract 

A completely decentralized adaptive fuzzy control 
procedure is proposed for a class of large-scale nonlinear 
systems with strong interconnections. The feedback and 
adaptive controller for each subsystem depend only upon 
local measurements to provide asymptotic tracking of a 
reference trajectory. The nonlinearity is approximated 
using a fuzzy logic method which transforms the 
nonlinear problems to linear ones. The unknown 
interconnections and external disturbances are predicted 
using the orthogonal basis of Hermite. The result of this 
prediction is taken into consideration in the synthesis of 
the control law. The proposed controller is combined 
with a Lyapunov analysis to ensure the stability and 
convergence. The proposed adaptive algorithm is tested 
on simulation example to control an interconnected 
double inverted pendulum. 

Keywords : Decentralized Control, Fuzzy Controller, 
Interconnected Nonlinear System, Orthogonal Basis. 

1. Introduction 

Decentralized adaptive control systems often arise from 
various complex situations where there exist physical 
limitations on information exchange among several 
systems for which there is insufficient capability to have 
a single control controller, and due to the physical 
configuration and high dimensionality of interconnected 
systems a centralized control is neither economically 
feasible nor even necessary. Therefore, the decentralized 
scheme is preferred in control design of large-scale 
interconnected system [1, 2, 6, 14]. To control a large- 
scale system, one essential problem is how to handle the 
interactions among different systems. Intensive research 
has been devoted to the observer design for large-scale 
systems. Uncertainties in a large-scale system require the 
adaptive decentralized technique, for which many 
decentralized adaptive schemes have been developed, 
including the model reference adaptive control [1,3], and 
nonlinear control with a special class of interconnections 
[5]. 

These approach focus on stabilization, where the 
dynamics of subsystems are assumed to be known or to 
be linear with a set of unknown parameters. However, in 



practice, large-scale systems may contain significant 
uncertainties, and/or with unknown parameters in 
nonlinear forms and unknown structures. 

Fuzzy logic control as one of the most useful approaches 
for utilizing expert knowledge has been an active field of 
research the past decade [8, 9, 10]. Fuzzy logic control is 
generally applicable to plants that are mathematically 
poorly modeled and where experienced operators are 
available for providing qualitative guidance. The most 
important advantage of Fuzzy-logic-control schemes lies 
in the fact that the developed controllers can deal with 
increasingly complex systems and controllers without 
precise knowledge of the model structure of the 
underlying system dynamic. Recently there have been 
significant research efforts on these issues in fuzzy 
control system [8, 10, 11] but these approaches work 
only for large-scale systems with known or linear 
dynamics with a set of unknown parameters and 
interconnection bounds. In practice, however, not all 
states are usually available. 

Based on fuzzy logic control and simple adaptation laws, 
this paper presents an approach which can easily tackle 
the output tracking control problem of a class of large- 
scale nonlinear system with unknown interconnections. 
The adaptive approach is provided which identifies the 
isolated subsystem dynamics to produce a stabilizing 
controller. This approach ensures asymptotic tracking 
using only local measurement. The interconnection terms 
and external disturbances are predicted online by using 
orthogonal basis projection. The stability and 

convergence are guaranteed using the Lyapunov 
approach. 

The paper is organized as follows: section 2 describes the 
problem under investigation; in section 3, we introduce 
the indirect approach and the stability analysis. 
Simulation results are given to demonstrate the 
effectiveness of the proposed approach in section 4. 

2. Problem Formulation 

Consider a class of nonlinear interconnected SISO 
subsystems S t (i = 1,2,..., N) described as follows: 
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\n. =fi(Xi ) +g i (x i )u i +A i (x) 

= x i,l 



( 1 ) 



where x = [X[ r ,X 2 ,...,x^] T ,x ( - e 9? is the global state 



vector. u t ( t ) <e 9? and y t ( t ) e 91 are respectively the 
input and output of the subsystem S t . The functions 
f t (.) and gf.) are unknown and nonlinear. 
Aj-Qc) e 91 represent interconnections among unknown 
subsystems ( i = 1 , 2 ,..., N ) and external disturbances. 

The tracking error for subsystem S i is defined by: 
e i ~ yir ~ yi • The goal is to design an adaptive control for 
each subsystem using only states which will cause the 
output y t to track a desired output trajectory y ir in 
presence of strong interconnections and disturbances 
using only local measurements. 

We have to make the following assumptions: 

A.l The vector x- is assumed to be available for 

measurement, this ensure the observability of each 
subsystem S t . 

A.2 The controllability is ensured if the functions 
^ 0 (i- 1 , 2 ,...,N ), for each subsystem 5^ . 
Subsystem (1) can be expressed as: 

Si = A/ofi + b i [ fi (*i ) + 8i (Xi K + A ; (*)] (2) 



where 





0 


1 


0 • 


•• 0 
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•• 0 
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0 • 
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> b i = 


0 
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0 • 


•• 0 










(d,-\) 

£i, 1 — yi ’ Gi, 2 — Tz ? ? — f z 



Assume that the given reference y /r is bounded and have 
up to d t - 1 bounded derivatives. The reference vector is 

denoted as: Y ir = [y ir , y ir , y ir ,..., y^ rV} ] T . Define the 
tracking error of the i th subsystem as: e t = . Then 

the error vector of the i th subsystem is given by: 

ii =[e { , e„ • 



It is desired that the output error of the I th subsystem 
follow: 

e \ dl) + k i,di-\ e i dl 1} + ••• + k i0 e i = 0 

here the coefficients are picked so that each: 

L = Si l +k i4 x gf l 1 +... + k 1 



i,0 



(4) 

has its roots in the open left-half complex plane. 

If each subsystem S t is known and taken independently 
without disturbances; (\(x\,x 2 ,...,x N ) = 0 ) , then the 



primary control should be designed to have the following 



idealized control law [13]: 

u * = — 7 —i~fi (h ) + Kf e t + yjy ) 

gife) 



(5) 



In practice the function (/-(x-) and g;(x- ) ) are 
unknown, nonlinear and there is any restriction 
conditions about the nonlinearity. The terms A ; (x)are 
unknown and composed of disturbances and 
interconnections between the subsystems. In this case, 
the determination of the primary control law is very 
difficult to implement. To overcome this problem, we 

propose to estimate the functions ( f t (x • ) and g;(x- ) ) 

by using the fuzzy logic technique. The unknown term 
A;(x) which regroups the interconnections, high order 
and disturbances is predicted using the orthogonal basis 
of Hermit. The result of this prediction is taken into 
account to synthesise the auxiliary signal. 



3. Adaptive fuzzy decentralized control 

The process of fuzzy logic approach includes three steps: 
fuzzification, fuzzy inference and defuzzification. The 
fuziffication consist to represent the crisp value into 
grades of membership for linguistic terms of fuzzy sets. 
The inference problem is to determine the value of the 
output based on the knowledge base and rule base. The 
defuzzification consists of converting fuzzified output 
into crisp value. 

Take an universal fuzzy system /^(x- / 6 f ) with 

X; €= U Xm for some compact set U x , to approximate the 
uncertain term //(A/) where 6 i contains the tunable 

parameters. Here, the linearly parameterized fuzzy model 
[ 8 ] is employed in the approximation procedure. Then we 
replace f t (x • ) and ^ (x^) by the approximated functions 

fMi/Q-n ) and sMiKln) respectively. 

The fuzzy system is described by a collection of fuzzy 
rule if-then as follow: 

: if x] is F t \ and ... if xf is F l in _ then (p\ is C\ 
where F- and C\ are fuzzy sets respectively in U i and 

v*. 

The approximated functions /• (x . 1 0_ iX ) and 

Si (*i I —i 2 ) can expressed as: 

fi(x i l6 a ) = e^ l (p i {x i ) ( 6 ) 

8i(Xi IO i2 ) = e^ il <p i {x i ) ,(i = 1,2,-, N ) 
where 9 ik = [d] k ,0f k ,—,0' k ]e 9 l m ‘ is a parameter vector 

(k = 1 , 2 ) and ^(x ; ) = \cp] (x,.),^, 2 (x l ),...,<p' n ' (x,)] e tR m; 

is a regressive vector. We have using the singleton 
fuzzifier, center average defuzzifier and product 

inference to calculate the regressor cp\ (x • ) : 



14 



Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 8, Issue 3, ICGST, Delaware, USA, January 2009 







z n 

l 7=1 



(7) 



where p pt (.) is the membership functions for 1 < l < nil 



( mi is the number of rules) and 1 < j <n t . 

In this paper, we propose the decentralized adaptive 
fuzzy controller defined as: 



1 



-(-.DixJd^+Kje, 



( 8 ) 



Uhl On) 

+ y\r +U is +U ih ) 

u is and u ih are the auxiliary controllers. u is is introduced 

to perform the main control action and u ih is synthetised 

in order to compensate the effect of interconnections. 
According to the universal approximation theorem [4], 

there exist optimal approximation parameters 0* x and 0* 2 
such that fiiXj/Ol |) and g t (v • / 0 [ 2 ) can respectively 
approximate f (v ■ ) and g^.) as best as possible. 

Define the optimal parameter vectors and fuzzy 
approximation error: 



9 n = argmiivpup t e9i „, |/,(x,.) - /,(x, > O n ) 



On = argminpup^,, |g,-(x ( .) - g,(x,. 1 0 i2 )\ 

O i2 &Cl i2 1 

where Q iX and Q i2 are the convex compact sets which 

contain parameter sets for 6* x and 0* i2 respectively: 

n a =W-tr(0n0l)<M a ) 

^■n = 'fa : tr (0n0 T n) ^ M i2 1 
define 

V i (x i ) = f i (x i )-fi(x i 'On) , 

A 8i(x i ) = g i (x i )-g i (x i /d- 2 ) 

which denote the minimum approximation error. 
Substituting (7) into (1), the tracking error dynamic 
equation can be written as: 

= K i e i + b t [/,- ( x ) - f t (x,. / 0 n ) + (g,- (x ; ) 

- 8i (A ; , 0 i2 )H - u is - u ih - A,- (x)] 
where 



(9) 



( 10 ) 



( 12 ) 
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•• 0 
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•• 0 
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•• 1 


> b ,= 
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-Kto 


Kt i 




•• K idi 





(13) 



from (6) and (1 1), equation (12) can be written as: 
<?,. = K i e i + b, [< b\(p i (x ( . ) + $> T n<Pi (x ; )u t 
+ A/,- (i) + &gi (A; )m t ~ u is -u ih - A; (x)] 
where the i th parameter error vector is defined as: 
®il = 0;\ - 0 n , 0 ;2 = 0 i2 - On 



(14) 



The problem of the interconnections is relaxed by 
introducing the orthogonal basis functions [11]. We 
project the interconnection terms on the Hermite basis 
functions, A; (v) becomes: 

(15) 



A; (x) = r* T S( (t) + W A . 
where T* is a q t vector such that : 

r* 2 ••• r*^] r and w A . is the projection error. 
S t (0 is the Hermite basis of dimension q t . 

Let T t be the estimate of T* , u is is chosen as: 
u is =-TjS({t) (16) 

From equation (15) and (16), the error dynamics can be 
represented as: 

k( = A ; e,. +b i [®J l <p i (x i ) + <!>n<p i (x i )u i 

+ A/) (a) + A gi(X()U( - u ih (17) 

~rfS(( t )- vv A ] 

where f, = r*-T,. 

Throughout this section, we need the following 
assumption: 

A.3 There exists a positive function 0 < W i (x ■ ) such that: 

T 

j|H|^(x ; .) Vl<i<JV (18) 

0 

with W( = max(|® ; |,|w A . |) and C 0 ( = A/ ; (x) + Ag i (x i )u i . 

The following update laws are defined for the 
decentralized indirect adaptive controller: 

(19) 

( 20 ) 

( 21 ) 



On = 7,i PA<Pi( A,) 

On =0n e J Pi b i<Pi(h) u i 

1 T , 

u ih = P, b , 

2r ; 



(11) ^i=PA, e J Pi b i S i(t) 



( 22 ) 

u is =-TjS({t) (23) 

where rj iX , q i2 and r/ A are fixed adaptive gains. 

Theorem: Consider the nonlinear subsystem (1), if there 
exists matrices p t = pj >0 satisfying the Riccati- 
equations: 

1 2 

A T i Pi + P,A, + Q, - P(b ( ( j)b[ p( = 0 

r i P 

where Q t : = Qj > 0 , p > 0 is a prescribed attenuation 
level, and r { is a positive parameter to be designed 
subject to 2 r { < p 1 , while the control (8) is adapted 

according to (19)-(23), then the following properties are 
guaranteed for i - 1 , 2 ,..., N : 

(i) All the variables of the closed-loop system 
are bounded 

(ii) The performance tracking is achieved. 
Proof: Take the error dynamic equation (14) and 
consider the Lyapunov function candidate 
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(24) 



1 / 1 . / j 

v t =-e< Pi e i + — ®ii®n 
2 Pa 

+ — 0 r 7 0 +— !— f- r f- 

^ ^ i2^i2 ^ ~ x i x z 

»7/2 2? 7a ; 

where each /?, e is a positive definite and 

symmetric matrix. 

The time derivative of v t along the error trajectory 

1 -t 1 t 2 : 

v i =~ e i Pi e i +~ e i Pi e i + — 



Pa 

2 T ’ 1 ~ T ~ 

+ — ®f 2 ®i 2 + — r f r, 

Pa Pa, 

Substituting (14) into (25) we obtain: 
Vi = Lf ( A I Pi + pAi) e i + 

e] PA [<&]\(Pi (x t ) + c t>J 2<Pi (x. )«,. 

2 T ■ 

- u ih + co i - w a, ] + — 4 J (l O,| 

Pa 



(25) 



(26) 



+ — ®; 2 ®/ 2 +— r/r,. 



1 -r- 

r rT 

Pa, 



2 

Pa 

Applying parameters adaptation laws (19)-(23), yields 

Vi ^ ~\ej Q,e, - \-e] Pi b,(-\- —)bf p,e, 

2 2 P r i 

+ ej Pib;W , - ej Pib iUih 

^-\ejQ iei +^-(eJ Pl b ,) 2 

2 2 r,- 



■ \{ej PA ) 2 + ej pAw - ef p,- *>,•«,* 

P 



(27) 



1 



1 



^ -4-ef 2i«i --V(«r Pibi) 2 +ej PibjW, 

2 p 

^ ~ej Q,e, - (—e? Pi b, - ypw ; ) 2 

2 p 2 

+^(pw,) 2 

Consider the composite system Lyapunov candidate 

N 

V = ^ where each s { > 0 , 

i=i 

v < <2,e, +L I p\v, 2 ] 



(28) 



1=1 

Denoting: 

W = diag(^ Wl ,...,Jj; WN ), 

Q = diag(e l Q l ,...,£ N Q N ) and E = (e u ...,e N ) (29) 

one has 



V <--E T QE + -p 2 W T W 
2 4 



(30) 



Integrating the above inequality from t - 0 to T yields 

T T 

V f ""’" 2 
2 



T 

2 V(T) + J \\Efdt < 2V (0) + -p : 



dt 



(31) 



from V(T) and assumptions, it is clear that 

E(t)eQ E =\E/E T QE<2V(0) + -p 2 max(w i )\so, all 
[ 2 i<i<v J 

the variables are bounded. 



4. Simulation results 

To show the performance of the presented controller, we 
consider the double inverted pendulum [7]. 

The equations which describe the motion of the 
pendulums are defined by: 



*11 =*12, 



{ m l gr kr 



kr 



*12 = H*- - — )sin(v n ) + —(l - b) 
J\ 4j x 2 h 

u 1 kr 2 . . 

+ — + — sm(-r 2 i) 

J\ 4 Ji 



X 2l ~ X 22 j 



.m^gr kr 2 . . kr 

*22 =(— — )sin(x 21 ) + — ( 1 -b) 

72 4 j 2 2y 2 

w 2 kr 2 . , 

H — 7- H — sin(x 12 ) 



(32) 



(33) 



72 4 72 



where and x 2l -(j) 2 are the angular 

displacements of the pendulums from vertical. The 
parameters m l = 2 kg and m 2 = 2.5kg are the pendulum 
end masses, jj = 0.5 kg and j 2 = 0.625 kg are the 
moments of inertia, the constant of connecting spring is 
k = 100 N/m, the pendulum height is r=0.5m, the natural 
length of the spring is l = 0.5m and the gravitational 
acceleration is g = 9. 81m/ s 2 . 

The distance between the pendulum hinges is defined as 
b - 0.4m (with b <1 in this example, so that the 
pendulum links repels each other when both are in the 
upright position. The control parameters for simulation 
are chosen as follows: 

rj n = 0.001, 77 /2 = 0.001, // A = 0.0001, Q t -diag{ 10,10) , 

A 2 

and each A- so that L(g) = g +4^ + 4 has roots at 
(-2, -2) . 

In simulations we are chosen the same conditions as in 
[12]. 

(*n,*i2 5*2i,*22) - (n / 3,0,—^ / 3,0) ), 0 iX = 0 i2 = 0 2 5 X i 

and p = 0.1. 

We define three fuzzy sets for component of each 
x x =(x u ,x l2 ) and x 2 =(x 2l ,x 22 ) with labels A 1 ^ ; , , 

, A^ and A^ characterized by: 



^A 1 . .(*iy) = ex P(-(*i/+ 0 - 8 ) ) 



^a 2 = ex P(-(*i/ + °- 4 ) ) 



0 



0 



^A 3 . (* 17 ) = ex P(-(*y) ) 
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/V ( X ij ) = eX PH*y - °' 4 ) 2 ) 

x u 

(X U ) = ex PH*y - °- 8 ) 2 ) 

X ij 

with n t j = 5 , j = 1,2 and i - 1,2 . 

Defining 25 fuzzy rules, in the following linguistic 
description: 

r\ 1) :if x-j is and x i2 is A J x 2 then y\ is C\ (34) 

25 2 

Denoting D t = Y FT ju z (*#. ) , 

a ^ 

/=i *=i 

(Pi O; ) = O, )/Dj , (Pi, 2 VO, , • ■ ^,,25 ( x i )/D; ] r 

In order to test the proposed approach, different 
simulations are studied. For each, the reference signal is 
taken equal to zero. 

For the first and second simulation the system is not 
affected by a disturbances and the functions 

fi (x • ) and g i (x i ) are unknown nonlinear. (x t ) and 

gi(x t ) are approximated by a fuzzy logic approach. We 

note that for the first simulation (figure. 1) the 
interconnections are supposed know. 

In the second simulation (figure. 2) the unknown 
interconnections are predicted using the orthogonal 
Hermite basis truncated after three terms. 

It can be seen from figure. 1 and figure. 2 that a good 
tracking is obtained which prove the effectiveness of the 
predictions. 




time (sec) 

Figure 1. The behaviour of the states: the interconnections are know 




time (sec) 

Figure 2. The behaviour of the states: the interconnections are 
predicted 



The third and fourth simulations are realised under the 
same conditions as for the first and second simulations 
respectively. A Gaussian noise (figure.5) is injected to 
the out puts of each subsystem. 

For the third simulation, the Gaussian noise is not taken 
into consideration, figure.3 show that the convergence is 
less good. When the Gaussian noise is added to unknown 
interconnections and predicted using the orthogonal 
basis; a good tracking is illustrated on figure.4. 





Figure 3. The behaviour of the states: the system is perturbed and the 
interconnections are know 




time (sec) 

Figure 4. The behaviour of the states: the system is perturbed and the 
interconnections are unknown 



20 
18 
16 
14 
12 
10 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

time (sec) 

Figure 5. Gaussian noise 




5. Conclusion 

In this paper, we have presented an adaptive output- 
feedback fuzzy decentralized control for a class of large- 
scale nonlinear systems. In the proposed method, fuzzy 
logic technique is used to estimate the part of the system. 
The interconnection terms and external disturbances are 
predicted by projecting them onto orthogonal basis. The 
adaptive control laws are derived from lyapunov analysis. 
A good convergence and tracking are obtained without 
any restrictive conditions on the interconnection terms 
and disturbances. 
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Abstract 

This paper describes the Load Frequency Control (LFC) 
of two area interconnected reheat thermal system using 
conventional Proportional - Integral (PI) controller and 
Fuzzy Logic Controller (FLC). The system is 
incorporated with governor dead band, generation rate 
constraint non-linearities and boiler dynamics. The 
conventional PI controller does not yield adequate 
control performance with the consideration of non- 
linearities and boiler dynamics. To overcome this 
drawback Fuzzy Logic Controller has been employed in 
the system. The aim of FLC is to restore the frequency 
and tie line power in very smooth way to its nominal 
value in the shortest possible time. Time domain 
simulations are used to study the performance of the 
power system. System performance is examined 
considering 1% step load perturbation in either area of 
the system. 

Keywords: Area Control Error, Fuzzy Logic Controller, 
Governor Dead Band, Load Frequency Control, 
Proportional-Integral. 

1. Introduction 

With an increasing demand for electric power, the 
electric power system becomes more and more 
complicated. Therefore the supply of electric power with 
stability and high reliability is required. Normally, the 
power system operates in normal state which is 
characterized by constant frequency and voltage profile 
with certain system reliability. For a successful operation 
of power system under abnormal conditions, mismatches 
have to be corrected via supplementary control [4]. 
Automatic Generation Control (AGC) or Load Frequency 
Control is a very important issue in power system 
operation and control for supplying sufficient and reliable 
electric power with good quality. An interconnected 
power system can be considered as being divided into 
control areas, which are connected by the tie lines. In 
each control area, all generators are assumed to form a 
coherent group. The power system is subjected to local 



variations of random magnitude and duration. For 
satisfactory operation of a power system the frequency 
should remain nearly constant. The frequency of a system 
depends on active power balance. As frequency is a 
common factor throughout the system, a change in active 
power demand at one point is reflected through out the 
system. Mostly the boiler system effects and the governor 
dead band effects are neglected in the load frequency 
control studies for simplicity. But for the realistic 
analysis of system performance, these should be 
incorporated as they have considerable effects on the 
amplitude and settling time of oscillations [9], [11]. 

Due to unnecessary error in the past, conventional PI 
controller does not provide adequate control performance 
with the consideration of non-linearities and boiler 
dynamics. The difficulty in obtaining the optimum 
settling time of previously said controller is mitigated by 
using FLC, which gives the opportunity to describe the 
control action in qualitative term and symbolic form. 
Literature survey shows that [8], [13-15] only a few 
investigations have been carried out using FLC to LFC 
with the consideration of non-linearities and boiler 
dynamics. It has been discussed that the implementation 
of FLC has greatly improved the performance of the 
controller “without negatively affecting the consumers’ 
quality of supply”. The remainder of the paper is 
organized as follows: Section (2) focuses on the transfer 
function modeling of two area thermal system 
considering boiler dynamics and non-linearities for load 
frequency control studies. Section (3) emphasizes on 
conventional proportional-integral (PI) control strategies. 
Section (4) discusses the proposed fuzzy logic controller. 
Finally section (5) and section (6) presents the simulation 
results and conclusion respectively. 

2. System Investigated 

The detailed block diagram modeling of two area thermal 
power system, for load frequency control, investigated in 
this study is shown in figure (1) [5], [6]. Here, area 1 and 
2 comprising reheat thermal system with governor dead 
band, generation rate constraint non-linearities and boiler 
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dynamics. The nominal parameters of the system are 
given in Appendix. Matlab version 7.3 has been used to 
obtain dynamic response such as frequency deviation in 
area 1 (AF^, area 2 (AF 2 ), tie line power deviation (AP tie ), 
and change in steam flow in area 1 (AS f i), area 2 (AS f2 ), 
for 1% step load perturbation in either area of the system. 

2.1 Governor Dead Band 

Governor Dead Band (GDB) is defined as the total 
magnitude of a sustained speed change within which 
there is no resulting change in valve position. The 
Backlash non-linearity tends to produce a continuous 
sinusoidal oscillation with a natural period of about 2s 
[1]. The speed governor dead band has significant effect 
on the dynamic performance of load frequency control 
system. Describing function approach is used to 
incorporate the governor dead band non-linearity [3]. The 
hysteresis type of non-linearities is expressed as, 

• 

y = F(x, x) rather than as y = F(x ) (1) 

To solve the non-linear problem, it is necessary to 
make the basic assumption that the variable x, appearing 
in the above equation is sufficiently close to a sinusoidal 
equation, that is, 

x^ASinco 0 t ( 2 ) 

where, A is amplitude of oscillation 

CO 0 is frequency of oscillations 

— o ^ 

As the variable function is complex and periodic 
function of time, it can be developed in a Fourier series 
as follows, [3] 

• tv • 

F(x , x) = F° + N x x H — - x+ (3) 



As the backlash nonlinearity is symmetrical about the 

origin, F° is zero. For the analysis in this paper, 
backlash of approximately 0.05% is chosen [1]. From the 
above equation, for simplification, neglect higher order 
terms, the Fourier co-efficients are derived as Ni=0.8 & 
N 2 =-0.2. By substituting the values in equation (3) the 
transfer function for GDB is expressed as follows, 

• 0.2 * 

F(x,x) = 0.8x x (4) 

71 

2.2 Generation Rate Constraint 

In practice, there exists a maximum limit on the rate of 
change in the generating power. For thermal system a 
generating rate limitation of 0.1 p.u MW per minute is 
considered, [10], [13] i.e. 

• 

A P g <0.1 p.u.MW / min = 0.0017 p.u.MW / s (5) 

2.3 Boiler Dynamics 

Figure (2) shows the model to represent the boiler 
dynamics [7]. Boiler is a device meant for producing 
steam under pressure. The model is basically for a drum 
type boiler. An oil/gas fired boiler system has been 
employed in this study, since such boilers respond to load 
demand changes more quickly than coal-fired units [1], 
[3]. Drum type boiler is otherwise known as recirculation 
boiler which relies on natural or forced circulation of 
drum liquid to absorb energy from the hot furnace walls, 
called water walls for generating steam. 

The boiler receives feed water which has been 
preheated in the economizer and provides saturated steam 
outflow. Recirculation boiler make use of a drum to 
separate steam flow from the recirculation water so that it 
can proceed to the super heater as a heatable vapour; 
hence recirculation boiler are referred to as drum type 
boiler[2]. The changes in generations are initiated by 
turbine control valves and the boiler controls respond 
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with necessary control action, changes in steam flow and 
changes in throttle pressure, the combustion rate and 
hence the boiler output. 




Figure 2: Boiler dynamics 



3. Conventional PI Controller 

The proportional-integral (PI) controller has received a 
great deal of attention in the process control areas. It is 
used as a feedback controller which drives the plant to be 
controlled with a weighted sum of the error and the 
integral of that value. 

u x --K p • ACE X - K t J ACE x dt (6) 

u 2 =-K p ■ ACE 2 - K i | ACE 2 dt (7) 

Where, K p and Ki are proportional and integral gains 
respectively, ACE is area control error which defines “a 
quantity reflecting the deficiency or excess of power 
within a control area” [17] and ul, u2 are controlled 
output of the respective areas. 

The relative simplicity of this controller is a successful 
approach towards zero steady state error in the frequency 
of the system. When the integral term is combined with 
proportional controller it accelerates the movement of the 
process towards set-point and eliminates the residual 
steady state error. The conventional control strategy for 
the LFC problem is to take the error as the control signal. 
For conventional PI controller the gain K p and Ki has 
been determined using Integral Square Error (ISE) 
criterion. The objective function used for this technique 
is 

t 

/ = | ( AFj 2 + AF 2 + AP t f e ) dt (8) 

0 

where, 

AF = change in frequency and 
AP tie = change in tie line power 

On the basis of performance index curve as shown in 
the figure (3) feedback gains should be determined to 
achieve the optimality of system performance. The main 
goal of FFC in interconnected power system is to protect 
the balance between generation and consumption. 
Because of the complexity and multi variable condition 



of the power system, a conventional model may not give 




Figure 3: Performance Index Curve 



Using the ISE technique for PI controller, the optimum 
value of K p and Ki are found to be 0.04 and 0.25 
respectively. However, it is learned from the responses 
with PI controller that, the conventional model responses 
are unstable. This inference is further reinforced by fuzzy 
logic controller which is useful in mitigating a wide 
range of control problems. 



4. Fuzzy Logic Controller 

The concept of fuzzy logic was developed by 
Eotfi.A.Zadeh in 1965 to address uncertainty and 
imprecision which widely exists in engineering problems 
[12], [15], [16]. His process approach emphasized 
modeling uncertainties that arise commonly in human 
thought processes. The design of FFC can be normally 
divided into three areas namely allocation of area of 
inputs, determination of rules and defuzzifying of output 
into a real value. In this paper the proposed fuzzy 
controller takes the input as ACE and ACE, which is 
given as, [15] 

ACE } =AF i B i +AP tie (9) 



The method of fuzzification has found increasing 
applications in power systems. The applications of fuzzy 
sets signify a major enhancement of power systems 
analysis by avoiding heuristic assumptions in practical 
cases. This is because fuzzy sets could be deployed 
properly to represent power system uncertainties. 

The parameter Bi may be optimized, but here, chosen 
as, [13] 




1 



+ — 

R, 



( 10 ) 



The Block diagram of Fuzzy Fogic Controller is 
shown in figure (4). Membership Function (MF) specifies 
the degree to which a given input belongs to a set. Here, 
seven membership function have been used to explore 
best dynamic responses, namely Negative Big (NB), 
Negative Medium (NM), Negative Small (NS), Zero 
(ZO), Positive Small (PS), Positive Medium (PM) and 
Positive Big (PB). 

Fuzzy rules are conditional statement that specifies the 
relationship among fuzzy variables. These rules help us 
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to describe the control action in quantitative terms and 
have been obtained by examining the output response to 
corresponding inputs to the fuzzy controller. Rules are 
given in Table I. The rules are interpreted as follows, 



If ACE is NB and ACE is NS then output is PM. 




Figure 4: Fuzzy Logic Controller 




Defuzzification, to obtain crisp value of FLC output is 
done by center of area method. 



Table I 

Fuzzy Control Rules 
ACE 





NB 


NM 


NS 


ZO 


PS 


PM 


PB 


NB 


PB 


PB 


PB 


PB 


PM 


PM 


PS 


NM 


PB 


PM 


PM 


PM 


PS 


PS 


PS 


NS 


PM 


PM 


PS 


PS 


PS 


PS 


ZO 


ZO 


NS 


NS 


NS 


ZO 


PS 


PS 


PS 


PS 


ZO 


NS 


NS 


NS 


NS 


NM 


NM 


PM 


NS 


NS 


NM 


NM 


NM 


NB 


NB 


PB 


NS 


NM 


NB 


NB 


NB 


NB 


NB 



5. Simulation Results 

Simulations were performed to a two area reheat thermal 
system considering non-linearities and boiler effects. The 
system is simulated for a step load disturbance of 1 % on 
area 2. Due to this, change in dynamic responses of the 
system has been observed as shown in figure 5(a)-(e). It 
is observed that the FLC exhibits relatively good 
performances having very smaller overshoot and 
transient frequency oscillations. The conventional PI 
controller does not yield adequate control performance. 
Simulation result concludes that FLC yields much 
improved control performance than the conventional PI 
controller. 
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6. Conclusion 

This study is an application of FLC to load frequency 
control in the power system. In this work, dynamic 
behaviour of the frequency of each area, tie-line power 
deviation and change in steam flow of each area in the 
power system with two area reheat thermal power system 
considering non-linearities and boiler dynamics under 1% 
step load disturbance in area 2, is studied. Integral square 
error technique has been used to obtain conventional PI 
controller gains. The conventional PI controller does not 
yield adequate control performance with the 
consideration of boiler effects and non linearities. Seven 
triangular membership function has been used to design 
FLC. Here, ACE is used as error signal, to control the 
power system. Finally the simulation results of two-area 
system with FLC have been compared with conventional 
PI controller. The simulation results conclude that FLC 
yields fast settling time with less number of oscillations 
which advocates the smooth settlement of the quality 
power supply. Also, the result shows that the FLC yields 
much improved control performance when compared to 
conventional PI controller. 
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8. Appendix 

(a) System data: [1] 

P ri =2000 MW 
T ti =0.3s 
T gi =0.08s 
K ri =0.5 
T ri =10s 

K pi =120 Hz/pu MW 
T pi =20 s 
T 12 =0.086 
Ri =2.4 Hz/pu MW 
f=60 Hz 

Bi =0.425 pu MW/Hz 
i = 1 & 2 
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(b) Boiler (oil fired) data: [1] 



Ki= 0.85 
K 2 = 0.095 
K 3 = 0.92 
C b = 200 
Td= 0 
Tf= 10 
K ib = 0.03 
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Abstract 

In this paper, we present a novel maximum power point 
tracking method for a photovoltaic system consisting of a 
photovoltaic panel with a power electronic converter; the 
whole is feeding a battery. This maximum power point 
depends on the temperature and irradiation conditions. A 
robust control using a PI regulator is used to track this 
maximum power point. The synthesis of this regulator 
has been achieved by using Bode method. For having a 
transfer function of the system, we have used a small 
signal modeling. Satisfactory theoretical and simulated 
results are presented. The simulation gives good results. 

Keywords: Photovoltaic panel, boost dc/dc converter, 

PI regulator, Maximum Power Point Tracking. 

NOMENCLATURE 

MPPT maximum power point tracking 

PV photovoltaic panel 

MPP maximum power point 

1. Introduction 

Due to the high cost of solar cells, it is necessary that the 
PV module operates at its maximum power point. The 
maximum power produced by a solar cell changes 
according to the solar radiation and temperature. A PV 
module is a nonlinear generator. The typical power- 
voltage characteristic of photovoltaic panel is shown in 
figure 3. The MPP (V m , P m ) is reached where 

dPj dV = 0, P = VI being the PV power. In order to 
optimize the ratio between the output power and the 
installation cost, photovoltaic systems have to draw 
continuously a maximum power from the modules. 
Maximum power point trackers, commonly known as 
MPP-trackers, are systems forcing PV modules to operate 
at their maximum power. Therefore, various solutions are 
presented in the literature. Many MPPT systems use a 
microcontroller or a personal computer for implementing 
sophisticated algorithms [1-4], or even neural networks 
[5]. These systems ensure very high performances. 
However they are typically very expensive and often 



need a separated, stable power supply for its operation; 
therefore they are only suitable for high power 
applications [6]. Another algorithm is based on the 
searching of operating point which makes 

dP/dV =0. Since the sign of dP/dV gives the 
direction of the search; it is possible to determine the 
maximum power operating by continuously detecting the 
PV power and voltage. In recent years, many MPPT 
applications based on this searching algorithm have been 
presented [7]. In [8, 9] an analog controller is proposed, 
where a boost dc/dc converter is handed for having 

dI b /dV equals zero (Figure 1). In this method, in order 

to reduce the complexity of the system, we have 
considered that the battery is equivalent to a constant 
voltage E on one hand. On the other hand, the converter 
is assumed to be ideal. 

In this paper, we pick up again the work proposed in [9], 
and we consider that the battery is a constant voltage E in 
series with a constant resistance R b [10-12]. 

A block diagram of the proposed system is shown in 
figure 1 . A boost dc/dc converter, considered to be ideal, 
is used to interface the PV output to the battery and to 
track the maximum power point of the PV module. Then 

the controller must keep (dP/dV) equals zero. That is 
possible with action on the duty cycle (X (0 < a < 1) 

according to the solar irradiation X and to the temperature 
T. Duty cycle is a signal produced by a PI regulator. For 
synthesizing this regulator, we have developed a transfer 
function for the system using the small signal model. 
Coefficients K p and IQ of the PI regulator are obtained by 
frequency synthesis. 




Figure 1 : Block diagram of proposed system 
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The paper is organized as follows: after the Theoretical 
study and simulation procedure (sections 2 and 3), the 
theoretical results are presented in section 4 and the 
simulation results are proposed in section 5. Conclusions 
are presented in section 6. And finally, the references are 
given in section 7. 

2. Theoretical study 

2.1. Optimal operating 

The equivalent circuit of a PV module is shown in figure 
2. The PV module considered in this paper is the SM55. 
It has 36 series connected mono-crystalline cells. 

The relationship between the panel output voltage V and 
current I is given by authors in [13, 14] as follows: 

V + RJ 
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( 1 ) 



( 2 ) 



(3) 



and 



I, V, I os , T, k, q, K I? I sc , X, Isoi, E go , Y(=P), T r , I or , R : 

R s are respectively, cell output current and voltage, cell 
reverse saturation current, cell temperature in Kelvin, 
degree Boltzmann’s constant (1.381e' 23 J/K), electronic 
charge (1.602e 19 C), short-circuit current at 25 °C and 
1000 W/m 2 , short-circuit current temperature coefficient 
at I sc (Ki =0.0004 A/K), solar irradiation in W/m 2 , light- 
generated current, band gap for silicon (=1.12 eV), 
ideality factor (=1.740), reference temperature (=298.18 
K), cell saturation current at T r , shunt resistance and 
series resistance. 

The output power of PV panel is P=VI , at optimal 
point, we have: 

— =/+ v — = 0 =>— =-— 

dV dV dV V 



>/ = (V-RJp a Aexp[A(v + R/)]+^- 



(4) 



Where: A = 

^TN cell 

cells in the module. 



and N ce ii is the number of series 





Figure 2: PV panel equivalent 
circuit 



Figure 3: PV panel P-V 
characteristic 



we have measured the shunt resistance R sh and 
determined the constants in these equations using the 
manufacturer ratings under standard test conditions of the 
PV panel. 

2.2. Boost converter study 

Value L of the inductor required to ensure the converter 
operating in the continuous conduction mode (figure 5) is 



calculated such that the peak inductor current at 
maximum input power does not exceed the power switch 
current rating [15]. Hence, L is calculated as: 

. (!-«»)«» 



V , 

L> om 



(5) 



fs\^Lm\ 

Where f s , a m , AI Lm , V om and I om are respectively, 
switching frequency, duty cycle at maximum converter 
input power, peak-to-peak ripple of the inductor current, 
maximum of dc component of the output voltage, and 
dc component of the output current at maximum output 
power. 

Taking into account that the ripple of the PV output 
current must be less than 2% of its mean value [15], the 
input capacitor value is given by: 



c> 






( 6 ) 



0.02(1 - a m y „ J , 

Where V inm is the PV module input voltage at the MPP. 
When the boost converter is used in PV applications, the 
input power, voltage and current change continuously 
with atmospheric conditions. Thus, the converter 
conduction mode could change since it depends on them. 
Also, the duty cycle a is changed continuously in order to 
track the maximum power point of the PV array. The 
choice of the converter switching frequency and the 
inductor value is a compromise between the converter 
efficiency, the cost, the power capability and the weight. 

2.3. Boost converter model 

If the chopping frequency is sufficiently higher than the 
system characteristic frequencies, we can replace the 
converter with an equivalent continuous model (figure 6). 
We will consider, for that, the mean values, over the 
chopping period, of the electric quantities. The transistor 
had been replaced by a voltage source whose value 
equals its mean voltage. At the same, the diode hade been 
replaced by a current source. Thus: 

V T = E( 1 -a)+ R b I b and I b = (l - a)l L (7) 

Where I L is the mean value of inductance current. 

We deduce from the continuous model equations: 

C^ = I-I L (8) 

at 



L d -E = v 

dt 



(1 -a)E-R b I b 



(9) 



bin 



( 10 ) 



At optimal operating point, we have: 

I m I Lm 

■V.=(\-a m )E + R b I 

Jbm ~ (l ~~ a m Lm 

Where I m , I Lm , I bm , V m , a m are respectively the optimal 
values of I, I L , I b , V and a. 

2.4. Small signal model and transfer function 

After having put the converter "mean" equivalent circuit 
in equations, we can determine the transfer function of 
the system open loop for a small variation around an 
optimal operating point. Thus, we write q = q m + A q, for 

each quantity q in the set ,(X,I L ,I ,I b ,I sol ,y] 
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defining the operating point. Where: y(t) = (t). So, 

dV 

the system can shown by a functional diagram like that 
used in [16] for synthesizing the regulators. The system 
can be presented as follows: 




Figure 4: Functional diagram of system. 





Figure 5: Boost converter waveforms (a) continuous conduction mode (b) discontinuous conduction mode 




I 





(b) 



Figure 6: (a) Boost dc/dc converter, (b) boost converter "mean" equivalent circuit 



Where: 



G al (s) = K al 

G a2 (s) = K a2 



1 

T b^s + 1 
s + s + 1 

CI 2 S + CL^S + 1 



( 11 ) 

( 12 ) 



With: 

K (K,G-K 2 lE + R t I u Jl + R,G m ) 

K,+K 2 R t (\-a m ) 



K t +K 2 R h \ 


[1- 


a ) 

m ) 


1 




[1-' 


«»)] 



And: 

b K^LC , b _ K 1 CR b (l-a m )+K 2 L 

2 K l+ K 2 R h (\-a m y ' K 1+ K 2 R b {l-a m ) 
n n GR b (l-a m ) + GL 

2 l + GR b (l-a m y 1 l + GR b (l-a m ) 



K =l , ARJRJ m -VJ . K _ A(R,I m -VJ G 
’ R dm G + Rfi m ) ’ 2 R d J\ + R s G m ) 

C _ G rn q _ 1 [ 1 

l + R,G m ' m R dm R sh ’ 



K al and K a2 are constants gain for a given temperature T 
and solar irradiation X. 

R dm and G m are the optimal dynamic resistance and 
conductance of PV module. 

The values of K ai , K a2 , K i? K 2 , b 2 , bi, a 2 , ai, G, G m and 
R dm can be easily determined if we have know the values 
of L, C, R b , E, I m , V m , a m and the parameter values of PV 
module. 

The transfer function of the open loop used to 
synthesizing a PI regulator is: 

G 0 (s) = G«(s)G„(s) = K M K„— r ± (13) 

U 2 S + CL^S + 1 

3. Simulation procedure 

We have built our model by using Simulink Matlab. The 
bloc used for simulations is given by figure 7. In PV 
module block, equations (1-3) are used; and in block 
(converter + battery), equations (7-9) are used. 




Controller 



„ 1 

/„Aexp[A(V„ +«,/.)]' 
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Figure 7: Bloc diagram for system 

The proposed controller circuit that forces the system to 
operate at its optimal operating point under variable 
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temperature and irradiation conditions is shown in figure 
8. Thus, for any temperature and solar irradiation level, 
the proposed controller circuit is obtained as follows: 

• On one hand, we multiply the PV output signal 
current I by the PV output signal voltage V. We 
obtain the PV output signal power P who is 
derived in order to obtain the (dP/dt) signal. 

• On the other hand, we derive the signal voltage 
V, who is inverted. Thus, the signal 1 /(dV/dt) 
is obtained. 

The product of 1 /(dV/dt) by (dP/dt) signals, gives the 

( dP / dV ) signal, who is compared to zero. 

The resulting difference signal (error signal) is used as an 
input signal of the PI regulator. This PI regulator is used 
to deliver a duty cycle signal to the dc/dc converter 

corresponding to the condition: _ q . 

dV 




Ground 



Figure 8: Controller circuit 



4. Theoretical Results 

4.1. PV module parameters 

Measure of R sh . The shunt resistance R sh of PV module 
(Figure 2) is measured. This measure (Figure 9) is 
performed in dark room (I so i = 0) in laboratory (no wind) 
by ammeter and voltmeter. We have applied an external 
negative voltage to the PV panel (I D = 0). 




Figure 9: Measure of R sh of PV panel. 

The measures of V mes and I mes must be done quickly (no 
overheating of photocell). The mean value of (V mes / I mes ) 
gives R sh =6500Q. R s is considered negligible in front of 
Rsh- 



Determination of the others PV module parameters. 

The manufacturer rated values of the SM55 PV module 
considered in this paper under standard test conditions 
(irradiation X = lkW/m 2 , A.M. = 1.5 solar spectrum and 
cell temperature T = 25 °C) are: open-circuit voltage V oc 



= 21.7 V, short-circuit current I sc = 3.45 A, maximum 
power current I m = 3.15 A, maximum power voltage V m 
= 17.4 V and maximum power P m = 55 W. 

Using these electrical characteristics, the constants I or? y 
(= P) and R s in equations (1) and (2), obtained by 
programming on Matlab Software are 4.842 pA, 1.740 
and 0.1 124 Q respectively. 

4.2. PI regulator and converter parameters. 

The PI controller gain and the integral time constant 
obtained by frequency synthesis using Bode method are 
respectively K P =0.01 and (=1/Kj) = 1.8 ms. From 
equation (5), the boost inductance can be chosen as is L = 
ImH. From equation (6), the input capacitance can be 
chosen as C = 4.7pF. The battery voltage and resistance 
utilized in this paper are respectively E = 24V and R b = 
0.65Q. Fifty kilohertz switching frequency is adopted. 

4.3. Optimal Values at MPP. 

For different values of irradiance X and temperature T, 
the computation of the theoretical optimum quantities V m , 
P m and a m of P, V and a are assembled in table 1 . 





V m (V) 


P m (W) 




X = 100 W/m 2 and T = 270. 1 8 K 


16.77 


5.270 


0.3069 


X = 1000 W/m2 and T = 270.18 K 


19.64 


62.73 


0.2468 


X = 1000 W/m2 and T = 320.18 K 


15.65 


48.61 


0.3984 


X = 100 W/m2 and T = 320.18 K 


12.30 


3.714 


0.4912 



Table 1. Quantities V m , P m and a m for different values of X and T. 

5. Simulation results 

In order to demonstrate the effectiveness of the proposed 
method, some simulations have been carried out. 
Simulations were achieved with Matlab Simulink. The 
simulations studies were made to illustrate the response 
of the proposed method to rapid temperature and solar 
irradiance change. For this purpose, the irradiance X and 
the temperature T, which are initially 100W/m 2 and 

270.18 K, are switched, at 0.05s and 0.15s, to 1000 W/m 2 
and 320.18 K respectively [Figures 10 and lObis (a), (b) 
and (c)] and vice versa [Figures 11 and 1 Ibis (a), (b) and 
(c)], i.e., the solar irradiance changes from 1000 W/m 2 to 
100 W/m 2 at 0.05s and the temperature changes from 

320.18 K to 270.18 K at 0.15s. 

[Figures 10(a), (b) and (c)] and [Figures 11(a), (b) and (c)] 
give the simulation results for values of K p and Tj of 
regulator PI obtained by Matlab programming. After 
adjust of Tj (=0.01ms) we have obtained the [Figures 
lObis (a), (b) and (c)] and [Figures 1 Ibis (a), (b) and (c)]. 
Figures 10(a), 10bis(a), 11(a) and 1 Ibis (a) give the 
variation of duty cycle, figures 10(b), 10bis(b), 11(b) and 
1 Ibis (b) give the variation of the output voltage of PV 
panel and the figures 10(c), 10bis(c), 11(c) and llbis(c) 
give the variation of instantaneous PV power for a step 
change of temperature and irradiance. It is shown that, 
the PI controller brings the system into the maximum 
power point after a some oscillations and the steady state 
is then reached. At the steady state, it is shown that the 
average values of voltage, instantaneous power of the PV 
and duty cycle of a boost DC/DC converter are very 
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close to theirs optimal values V m , P m and a m obtained by 
Matlab programming. So, the relative errors for different 
average values varied between 0.05% and 0.13%. 

The simulations of the MPPT show that the system is 
stable. The oscillations about the computed optimal 



operating point are due to the switching action of the 
DC/DC converter. The transients between operating 
points are natural for a dynamic system which is 
controlled by a PI type controller. 



simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




Time (s) 



— simulated Theoretical 




Figure 10: Variation before adjust of Ti, of (a) duty cycle, (b) PV 
output voltage, and (c) instantaneous PV power for a step change 
on irradiation X and temperature T from 100 W/m 2 to 1000 W/m 2 
and 270.18 K to 320 K respectively. 



Figure 10 bis: Variation after adjust of Ti, of (a) duty cycle, (b) PV 
output voltage, and (c) instantaneous PV power for a step change on 
irradiation X and temperature T from 100 W/m 2 to 1000 W/m 2 and 
270.18 K to 320 K respectively. 
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simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




simulated Theoretical 




Figure 11: Variation before adjust of Ti, of (a) duty cycle, (b) PV output 
voltage, and (c) instantaneous PV power for a step change on irradiation 
X and temperature T from 1000 W/m 2 to 100 W/m 2 and 320.18 K to 270 
K respectively 



— simulated Theoretical 




Figure 1 Ibis: Variation after adjust of Ti, of (a) duty cycle, (b) PV 
output voltage, and (c) instantaneous PV power for a step change on 
irradiation X and temperature T from 1000 W/m 2 to 100 W/m 2 and 
320.18 K to 270 K respectively. 



6. Conclusion 

The output power delivered by a PV module can be 
maximized using MPPT control system. It consists of a 
power conditioner to interface the PV output to the load, 
and a control unit, which drives the power conditioner for 
extracting the maximum power from a PV array. In this 
paper, a MPPT system has been proposed and tested by 
simulations in Matlab Software. It follows the irradiance 
and the temperature level change rapidly and tracks the 
MPPT. The PI regulator used to control the boost DC/DC 
converter is synthesized by frequency synthesis using 



Bode method. For that, we have developed a transfer 
function of global model using a small signal method. 
Simulations show that the regulation is robust against 
disturbances. 

For a given temperature T and solar irradiance X, we have 
calculated the corresponding optimal values a m , V m and 
P m . The optimal values obtained by programming 
coincide with ones obtained by simulation. The 
simulation gives good results. 

Our proposed MPPT’s charge controller is easier and 
cheaper to build. 
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Abstract 

Acoustic and vibration analysis are well proven non- 
destructive techniques for the diagnosis of machine 
health. Control valves are key to quality plant 
performance. Hence, their life span is of utmost 
concern for the operators. Condition monitoring for 
control valves has already been widely researched 
by several authors using the aerodynamic noise 
prediction method. This paper describes a method 
to predict the condition of the control valve, using 
acoustic signal. Acoustic analysis has been performed 
by evaluating the dBA levels for determining the 
faults in control valve. Sub-band approach is used 
for the acoustic signals to analyze the properties 
such as amplitude, frequency and energy change. A 
relationship between the acoustic signature with the 
faults has been established. 

Keywords: Globe valve , Acoustics, Signal pro- 
cessing, QMF. 



1 Introduction 

Control valves play a vital role in the control and 
optimal performance of a myriad of industrial flow 
processes. Their precise maintenance of pressures, 
flow, temperatures, and levels in thousands of 
installations is a key to overall successful plant 
performance. However, such stringent requirements 
and operating conditions induce several faults and 
failures in control valves. For the past three decades, 
extensive research is being carried out on predicting 
failures in control valves. There are many ways for 
establishing such failure diagnosis systems. Compu- 
tational Fluid Dynamics (CFD) Modeling has been 
used to predict the globe control valve performance 
characteristics. The parameters of interest were flow 
coefficient and inherent flow characteristics [1, 2], 
Erdal and Anderson have used a commercial CFD 
toolbox to evaluate various numerical effects in 



calculating the complex flow through a geometrically 
simple orifice. In this both the pressure drop and 
the downstream flow were evaluated and compared 
with experimental data [3]. From the CFD analysis 
of steam valve, it was observed that in the partial 
opening condition a high pressure region revolving 
circumferentially was the cause of vibrations [4]. 

Importance of Fault Detction and Diagnosis has 
been mentioned by A. Asokan and D. Sivakumar in 
the paper ‘Fault Detection and Diagnosis for a Three- 
tank system using Structured Residual Approach’ [5]. 
The study on Performance Evaluation of Integrated 
Fault-Tolerant Technique has been carried out by 
Abulnaja and Saadi[6]. Modeling and Diagnostics of 
Inductions Machines with neural networks has been 
proposed by Tarek Aroui et. al. [7] . 

Online control valve diagnosis has also been widely 
accepted which uses digital positioners. Kaseda et. 
al. [8] have developed two algorithms for detecting 
failures in operating parts of a control valve and one 
for detecting faulty valve operation due to stick slip 
while Schruers and Bednar [9] used an expert system 
to predict faults occurring in the field. For training 
the experts system, the actual maintenance records 
of the valve have to be catalogued and stored. 

Hyashi et. al. [10] have given detailed mathemati- 
cal analysis of how the vibrations are generated inside 
the control valves. With regards to control valve 
noise, ISA technical committee has drafted a standard 
for hydrodynamic and aerodynamic noise prediction 
in control valves. These documents provide the 
experimental setup and the tools needed to evaluate 
the noise levels inside the control valve in detail. 
IEC standard 534-8 7, Control Valve Aerodynamic 
Noise Prediction, involves first estimating the sound 
power generated in the fluid inside the valve and 
piping due to the throttling process [11]. Then, the 
transmission loss due to the piping is subtracted 
to determine the sound level at a predetermined 
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location outside the piping. A computer model has 
been built to predict the noise in control valve with 
the help of empirical data collected from experiments 
at a cogeneration plant [12]. One important issue 
is the piping vibration generated due to the control 
valve. Miller looked at the main causes of the flow 
induced vibration of the piping system. These causes 
are usually the two-phase flow situations involving 
cavitation or flashing, standing waves within the 
pipe, vortex shedding and a high kinetic energy 
exiting an upstream control valve trim [13]. Piping 
vibration can be used to measure the aerodynamic 
noise generated by the control valve. By measuring 
the vibration generated one meter downstream of a 
valve, the sound pressure levels can be calculated [14]. 
Research has been carried out for investigating the 
noise environments inside butterfly valves [15]. Work 
also has been done to diagnose the problems related 
with motorized operated control valves. Sharif et 
al. have used a pressure sensor along with a valve 
diagnostic system developed by Fisher Rosemount 
to perform fault diagnosis. Control valve parameters 
such as valve travel, frictional force, available seat 
load and spring rate were also studied[16]. 

The aim of this work is to predict faults occurring 
due to mechanical vibrations in control valves, using 
acoustic approach. For this work, acoustic data is 
captured using a capacitive microphone. A sound 
proof enclosure is used to eliminate the noisy process 
environment. The experimentation is carried on two 
different set-ups of different line sizes, with two globe 
valves. 

The remainder of the paper is organized as follows: 
Section (2) focuses on the experiment details for gath- 
ering data from the test rigs. The experimental results 
are provided in the Section (3). Section (4) gives the 
findings obtained from the experimental results. 

2 Experiment Details 

To perform diagnosis on control valves, the require- 
ment is to meet the demand for relevant data. The 
prime aim is to identify the faulty signatures of the 
control valves from the acoustic data so collected. 
The important blocks of the experimental set-up as 
specified in the international standards are given in 
the Figure 1. 



For the acoustic signal capturing, the sensor used 
is a microphone which converts the sound energy into 
an electrical output. The output can be directly fed to 
the sound card of a personal computer and stored for 
future analysis. A sound proof enclosure is provided 
around the control valve so as to avoid reflections from 
the walls and surrounding obstructions. Initially a 







I. High Pressure Source 
3. Upstream silencer 

5. Temperature Measuring Device 
7. Test Specimen 
9. Test Chamber 

I I. Downstream Throttling Valve 
13. Microphone 



2. Upstream Throttling Valve 
4. Flow Measuring Device 
6. Pressure Measuring Device 
8. Test Section 
10. Downstream silencer 
12. Low Pressure Receiver 



Figure 1: Test Rig as per IEC standard 



Sound Level Meter (SLM) is used for validation pur- 
pose. IEC standards for control valve noise define the 
use of anechoic chambers for aerodynamic noise test- 
ing of control valves. As interest is more in capturing 
the noise initiated in the moving parts of a control 
valve, a sound proof enclosure is used. The micro- 
phone is kept as close as possible to the valve to cap- 
ture the noise generated by the high turbulence region 
at the vena contract a and the valve trim. 

2.1 Microphone 

An externally polarized free-field condenser micro- 
phone is used to measure the acoustic fields near the 
control valve. The microphone is placed perpendicular 
to the control valve axis. Figure 2 shows the locations 
of the microphone on the control valve. The disturb- 
ing effects due to the presence of the microphone in 
the field are minimal. 




Figure 2: Microphone Mounting 



2.2 Sound Proof Enclosure 

Experiments suggest that due to the diffractions and 
reflections from the nearby surfaces, the microphone 
shows higher pressure levels. Standards as discussed 
earlier require the use of anechoic chamber for mea- 
surement of sound pressure level. Since the focus is on 
mechanical vibrations along the stem and trim assem- 
bly, sound proof enclosure is used as shown in Figure 
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3. This reduces the effects of reflections. The en- 




Figure 3: Inner view of the Sound Proof Enclosure 



pneumatic actuation. Table 2 gives the description of 
the faults introduced in the valves during experimen- 
tation. 



Table 2: Control Valve Faults 



Valve 1 


Valve 2 


Increased spring stiffness 
Decreased spring stiffness 
Loose packing 
Hunting 


Loose Hand Wheel 
Loose Brass nut 
Loose Body nut 



closure is covered with sound proof material from in- 
side to prevent external noise from affecting the signal 
strength. Sufficient gaps are made in the assembly to 
allow the piping through the enclosure without affect- 
ing the measurement. The details are shown in Figure 
4. 



3 Results 

The results are presented for the experiments con- 
ducted on control valves mounted on two different 
set-ups. 




Figure 4: Complete Test Assembly 



2.3 Test Assembly 

The tests are carried out on two different set-ups. The 
pump specifications for both the set-ups are identical. 
However, for the first set- up the line size is 25 NB 
and for the other it is 50 NB. 

2.4 Control Valves 

Specifications of the globe valves used for the experi- 
mentation are listed in Table 1. 

Valve 2 is provided with a hand wheel which reduces 



3.1 Valve 1 on Set-up 1 

3.1.1 Valve characteristics 

Initially the valve characteristics are plotted. The 
installed valve characteristics for healthy and faulty 
conditions is as shown in Figure 5. 




Figure 5: Characteristics of Valve 1 



The differential pressure across the Valve with 
opening of the valve is shown in 6. 




Table 1: Valve Specifications 



Properties 


Valve 1 


Valve 2 


Valve Size 


25 NB 


25 NB 


Valve Coefficient (Cv) 


7.5 


4 


Actuator 


Pneumatic 


Hand Wheel 


Characteristics 


Quick opening 


Parabolic 



the stem instability observed in the valve 1 due to 



Figure 6: Differential Pressure across the valve 1 
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3.1.2 Frequency Spectrum 

Frequency spectrum plotted from the unprocessed 
acoustic signal is as shown in Figure 7. 

3.1.3 Sound Pressure Level 

Table 3 lists the peak dBA values from the recorded 
signal for valve 1 on the first set-up. The signals are 
recorded for the healthy as well as for the commonly 
occurring faults. The nomenclature for various 
conditions is as follows: 

H = Healthy, FI = Increased stiffness of spring, 

F2 = Decreased stiffness of spring, 

F3 = Loose packing, F4 = Hunting 



Table 3: dBA values for valve 1 on set-up 1 



Valve opening(%) 


H 


FI 


F2 


F3 


F4 


0 


79 


83 


88 


84 


105 


25 


80 


86 


84 


82 


98 


50 


79 


85 


84 


84 


- 


75 


77 


84 


83 


83 


- 


100 


76 


84 


83 


82 


- 



These faults are true replica of the commonly oc- 
curring faults in valves installed in many process in- 
dustries. There are three types of faults viz. exter- 
nally correctable, non-correctable or correction needs 
to be done by dismantling the valve and loop related 
faults. Representative faults are considered while se- 
lecting the faults. Hunting is loop related whereas 
stiffness change can be corrected externally. For loose 
packing correction can be done by dismantling the 
valve. 

3.1.4 Energy Distribution 

The percentile energy in each frequency band is given 
in Table 4. 



Table 4: % Energy for Valve 1 on set-up 1 



Frequency Band(Hz) 


H 


FI 


F2 


F3 


215-315 


2.75 


3 


3 


5 


430-540 


0.26 


2 


1.6 


4 


860-900 


1.27 


15 


16 


9.5 


1075-1200 


0.1 


3 


4 


3 




(a) Healthy 



Spectrum of signal 
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(c) Decreased stiffness 
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3.2 Valve 1 on Set-up2 

The result of the experiments carried out on valve 1 
when installed in a flow loop set-up 2 are as follows. 
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Figure 7: Frequency spectrum of valve 1 on set-upl 
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3.2.1 Frequency Spectrum 

Spectrum of the recorded acoustic signals is given in 
Figure 8. 



3.2.2 Sound Pressure Level 

Table 5 lists the peak dBA values for the valve 1 on 
set-up2 for healthy as well as faulty conditions. 



Table 5: dBA values for valve 1 on Set-up2 



Valve opening(in %) 


H 


FI 


F2 


F3 


0 


87 


87 


85 


85 


25 


86.5 


87 


84 


87 


62.5 


88.5 


87.5 


90.5 


92 


75 


91 


87.5 


89 


85.5 


87.5 


93 


90.5 


85.5 


86 


100 


87.5 


89.5 


84 


85 



3.2.3 Energy Distribution 

The percentile energy in each frequency band is given 
in Table 6. 




(a) Healthy 




Frequency 



Table 6: %Energy for Valve 1 on Set-up 2 



Frequency Band(Hz) 


H 


FI 


F2 


F3 


40-45 


6.20 


2.0 


4.50 


12.0 


85-90 


9.50 


5.0 


7.50 


17.5 


185-200 


8.00 


0.3 


2.00 


1.50 


500-525 


7.75 


0.1 


1.75 


1.50 



3.3 Valve 2 on Set-upl 

For the valve 2, the microphone is placed near the 
stem that is perpendicular to the direction of acoustic 
source. 



(b) Increased stiffness 




Frequency 



(c) Decreased stiffness 



3.3.1 Valve characteristics 

The valve characteristics for valve 2 for the healthy 
and faulty conditions are as show in Figure 9. 

The differential pressure across the Valve with 
opening of the valve is shown in 10. 



3.3.2 Frequency Spectrum 

The fourier spectrum of the unprocessed signal and 
the corresponding faults is shown Figure 11. The 
nomenclature for various conditions is as follows: 




Frequency 



(d) Loose packing 

Figure 8: Frequency spectrum of valve 1 on Set-up2 
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% Opening 



[—♦—Healthy Faultl -*-Fault2 Fault3 1 



Figure 9: Characteristics of Valve 2 




Figure 10: Differential Pressure across the valve 2 

H = Healthy, FI = Loose Hand Wheel, 

F2 = Loose Brass Nut, F3 = Loose Body Nut 
As the second valve is hand operated, the faults re- 
lated to valve body are mainly considered. These 
faults are observed in the environment where struc- 
tural vibrations are predominant. 

3.3.3 Sound Pressure Level 

Table 7 gives the details of the dBA values of the 
valve for healthy and faulty conditions. 



Table 7: dBA values of valve 2 on Set-up 1 



Valve opening(in %) 


H 


FI 


F2 


F3 


0 


69.5 


65 


64.5 


71 


22.2 


81.5 


66 


71 


70 


44.4 


89 


85 


82.5 


89 


66.6 


84 


75.5 


80 


75 


88.8 


84 


85 


87 


83.5 


100 


91 


85 


90.5 


90 



3.3.4 Energy Distribution 

The energy levels in different frequency bands are 
given in Table 8. 



xio 3 Spectrum of signal 

2.0 | ! 1 1 ! 



1.5 




0 1000 2000 3000 4000 5000 6000 7000 

Frequency 



(a) Healthy 




Frequency 



(b) Loose Handwheel 



Spectrum of signal 




Frequency 



(d) Loose Body Nut 

Figure 11: Frequency spectrum of valve 2 on Set-up 1 
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Table 8: %Energy for Valve 2 on Set-up 1 



Spectrum of signal 



Frequency Band(Hz) 


H 


FI 


F2 


F3 


215-430 


3.7 


4.1 


1.6 


2.6 


1000-1215 


3 


5.7 


11.1 


1.3 


2150-2365 


1.5 


3.2 


6.63 


1.6 


3225-3440 


2.4 


1.25 


0.83 


2.39 



3.4 Valve 2 on Set-up2 

The result of the experiments carried out on valve 2 
when installed in a flow loop set-up 2 are as follows. 

3.4.1 Frequency Spectrum 

Spectrum of the recorded acoustic signals is given in 
Figure 12. 



3.4.2 Sound Pressure Level 

Table 9 lists the peak dBA values for the valve 2 on 
set-up 2 for healthy as well as faulty conditions. 



Table 9: dBA values for valve 2 on Set-up 2 



Valve opening(%) 


H 


FI 


F2 


F3 


0 


60.7 


62.75 


62.4 


65.5 


25 


62.5 


61.75 


61.8 


68.7 


62.5 


86 


61.8 


64.9 


70.7 


75 


65.8 


65.5 


65.8 


71.5 


87.5 


63.25 


65.7 


71.3 


77 


100 


79 


80.3 


78.8 


82.9 



3.4.3 Energy Distribution 

The percentile energy in each frequency band is given 
in Table 10. 



Table 10: %Energy for Valve 2 on Set-up 2 



Frequency Band(Hz) 


H 


FI 


F2 


F3 


40-45 


13.5 


5.5 


4.5 


9.5 


85-90 


0.95 


2.0 


2.75 


0.35 


185-200 


0.70 


1.5 


1.25 


0.75 


500-525 


0.5 


0.3 


0.1 


0.75 



4 Conclusion 

The purpose of this study was to identify the mechan- 
ical faults in globe valves using the acoustic approach. 



m 2.0 

p 

i 

i 




Frequency 



(a) Healthy 




(b) Loose Hand Wheel 




Frequency 



(c) Loose Brass Nut 




(d) Loose Body Nut 

Figure 12: Frequency spectrum of valve 2 on Set-up2 
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Experiments with different aspects for fault identifi- 
cation were carried out. The acoustic data analysis 
showed a change in the energy and sound pressure 
levels as regards to healthy and faulty conditions. As 
observed, each fault generated a pattern different from 
healthy as well as other faulty conditions. Hence a 
pattern recognition technique can be used to identify 
the health status of the valve. The limitation of this 
technique is the requirement of a sound proof cham- 
ber. 
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Abstract 

This paper presents a space vector based Generalized 
Discontinuous Pulsewidth Modulation (GDPWM) 
algorithm applied to direct torque control (DTC) of 
induction motor drive without angle estimation using the 
concept of imaginary switching times. The proposed 
GDPWM based DTC retains all the advantages of 
conventional DTC (CDTC) moreover eliminates the 
complexity involved in conventional space vector 
Pulsewidth modulation (CSVPWM) algorithm. The 
proposed method results in reduced current ripple over 
CSVPWM based DTC at higher modulation indices and 
low switching losses at all modulation indices. To 
validate the proposed method, simulation is done in 
MATLAB/SIMULINK environment and the results of 
some popular discontinuous PWM (DPWM) based DTC 
are presented. Also, comparison of the total harmonic 
distortion (THD) in line current between the CDTC and 
different DPWM based DTC is made. The proposed 
algorithm over and above enhancing the performance at 
high modulation indices, switching losses of the inverter 
and computational burden involved in conventional 
GDPWM algorithm are reduced significantly. 

Keywords: DTC, CSVPWM, DPWM, GDPWM Imaginary 
switching times. 



1. Nomenclature 

R s , R r stator and rotor resistances 



L s , L r , L m stator and rotor self inductances, mutual 
inductance 

P number of poles 

v ds ,v ds d, q axes stator voltages 



^ ds ’ ^ ds 



d, q axes stator currents 



V'ds’V'ds 

co. 



d, q axes stator flux linkages 
rotor electrical speed in radians 



CQ sl slip speed 

CO e synchronous speed in electrical radians 

T e electromagnetic torque 

J inertia constant of the induction motor 

2. Introduction 

Research interest in high performance control strategies 
for induction motor drives has grown significantly over 
the last three decades due to some of their advantages 
like less maintenance, simple construction and 
mechanical robustness. Many researchers paid attention 
to find out different solutions for the induction motor 
control having the features of precise and quick torque 
response, and to reduce the complexity involved in FOC 
[1]. The CDTC has been recognized as viable solution to 
achieve these requirements [2-6] .There is no need of 
current controllers in CDTC, as it selects the voltage 
space vectors according to the errors in flux and torque. It 
includes flux and torque estimators, flux and torque 
controllers and a switching table. The switching state of 
the inverter is updated in every sampling period. With in 
this sampling period, the inverter does not change its 
state until the outputs of the controllers change. For this 
reason the switching frequency is usually not constant; it 
changes with the rotor speed, load and bandwidth of the 
two controllers and consequently steady state ripples in 
torque and flux are high. To solve these problems some 
research was done on this. Recently several PWM 
techniques has been developed to reduce the steady state 
ripple and to get constant switching frequency. One of 
such techniques is the CSVPWM [7-9]. The steady state 
performance of the drive with this method is better than 
the CDTC and also the switching frequency remains 
fixed at a constant value. Though CSVPWM based DTC 
could solve most of the problems in CDTC it still has 
some drawbacks significantly at high modulation indices. 
Also, the CSVPWM suffers from the drawback of 
computational burden and increase in memory 
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requirement. In this paper an attempt is made to avoid the 
cumbersome calculations thereby eliminating the need of 
reference voltage vector position and sector calculations 
using the concept of imaginary switching times [10] -[12] . 
To improve the performance at high modulation indices, 
in recent years several DPWM are reported in [14] -[17]. 
Two popular approaches to PWM generation are the 
triangular comparison approach and the space vector 
approach. Performance of a particular DPWM method is 
modulation index dependent and so an attempt was made 
in [14], [15] to generalize the developed DPWM methods 
which lead to the development of generalized DPWM 
(GDPWM) algorithm based on the triangular comparison 
method [14]. 

The main objective of this paper is to propose a space 
vector based GDPWM algorithm for DTC induction 
motor drive using the concept of imaginary switching 
times, which yields high performance drive with reduced 
complexity. Section 2 presents a survey on the related 
work. Section 3 gives brief introduction to the CSVPWM 
algorithm. Section 4 introduces CDTC. Section 5 
explains the proposed technique. Section 6 presents the 
simulation results of different cases of the proposed 
technique. Finally Section 7 summarizes the present 
research work stating its merits and demerits. 

3. CSVPWM Algorithm 

With a 2-level, 3-phase voltage source inverter (VSI) 
there are eight possible switching states, among which 
six are active vectors and two zero vectors. It can be 
shown that all the six active vectors can be represented 
by space vectors given by (1) forming a regular hexagon 
and dividing the space plane into six equal sectors 
denoted as I, II ,111, IV, V, VI as shown in Figure 1. 

2 J(k-l£ 

V k =~V dc *e K k =l,2,...,6. (1) 



v 3 (0 1 0) V 2 (l 1 0) 




Figure 1. Switching states and corresponding voltage vectors of a three 
phase converter. I, II, III, IV, V, VI are the sectors. 

In the CSVPWM strategy, the desired reference voltage 
vector is generated by time averaging the two nearby 
active voltage states and two zero states in every sub- 
cycle^ . For a given reference voltage V REF at an an gle 
a with reference to V j in first sector, the volt-time 
balance is maintained by applying the active state Vj , 
active state V 2 and two zero states V 0 and V 7 together 



for durations, Tj,T 2 m\d ^respectively, as given in (2) 
[3], [10], [14]. 

T t = M * S 1,1(60 ~ a) * T s (2.1) 

Sin60° 

T 2 =M* — — —— * T s (2.2) 

Sin60° 

T Z =T s -Tj-T 2 (2.3) 

Where ‘M’ is the modulation index, given by 

m=2^. 

2V dc 

Also in CSVPWM strategy, the total zero voltage vector 
time, T z is divided equally by the two zero vectors 

V 0 and V 7 . Also in this method, the zero voltage vector 
time is distributed equally at the start and end of the 
sampling period in a symmetrical way. Thus, CSVPWM 
uses V 0 ,V k ,V k+1 ,V 7 ,V k+1 ,V k ,V 0 sequence for odd 

sectors and V 0 ,Vk+l ^7 >Vk+l >^0 sequence foe 

even sectors respectively. 

4. Principle of CDTC 

The principle of CDTC is that rapid instantaneous torque 
and flux control can be achieved quickly by changing the 
magnitude and position of stator flux linkage space 
vector. Decoupled control of torque and flux is achieved 
by controlling the orthogonal components of stator flux 
space vector. Neglecting the stator ohmic drop this can 
be made by properly selecting the inverter voltage states. 
An optimum switching table is constructed for selecting 
the appropriate voltage vectors to increase/decrease 
torque and flux. At every sampling period, based on the 
error of hysterisis torque and flux controllers an 
appropriate voltage vector will be selected so as to 
maintain the torque and flux error within the hysterisis 
band. 

In CDTC since only one vector is applied for the entire 
sampling period, so for small errors the motor torque may 
exceed the torque limit. Also because of the use of 
hysterisis controllers significant ripple in torque and flux 
is observed in steady state due to bang-bang behavior. 

5. Proposed GDPWM Based DTC 

The proposed space vector based GDPWM algorithm 
based DTC induction motor drive combines the 
principles of CSVPWM algorithm, concept of imaginary 
switching times, GDPWM algorithm and DTC. 

Concept of Imaginary Switching Times: The concept of 
applying imaginary switching times greatly reduces the 
complexity involved in CSVPWM algorithm in terms of 
reference vector position and sector calculation and 
memory required. With this approach the procedure of 
identifying the position of the reference voltage vector 
and sector estimation are eliminated. Assuming that the 
reference voltage vector is lying in first sector, according 
to the CSVPWM algorithm it can be best constructed in 
an average sense by the two of its nearest active voltage 
vectors (yj f V 2 ) and two zero voltage vectors (V^V/) in 
the specified sequence for a specified time Tj,T 2 ,T 0 ,T 7 
as in (2) respective [4]. It can be observed that actual 
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switching times Tj,T 2 are a dependent and so to 
eliminate the reference voltage vector position 
dependency the concept of imaginary switching times is 
applied to CSVPWM algorithm [9-13]. 

The active vector switching times Tj,T 2 are defined in 
terms of the imaginary switching times as given in (3) 

T l =T an~ T bn ( 3 -l) 

T 2= T bn~ T cn O- 2 ) 

The imaginary switching times proportional to the 
instantaneous values of the reference phase voltages are 
defined as 



T = 

± an 



T bn ~ 



\Y dc 



T = 



v 



bn 



dc 

\ 



\V dc J 



(4.1) 



(4.2) 



(4.3) 



Where, T s is the sampling period and V dc is the dc link 

voltage. For every sampling time, the maximum and 
minimum values of imaginary switching times are 
calculated and the active vector switching times are 
expressed in terms of the maximum and minimum 
switching times. 

The effective time is the duration in which the reference 
voltage vector lies in the corresponding active states, and 
is the difference between the maximum and minimum 
switching times. The zero voltage vector switching times 
is calculated from the switching time duration and 
effective time as, 



T Z= T s~ T eff where - 
T eff = T 1+ T 2 

After obtaining the switching times for the active states 
and zero states the reference voltage vector is synthesized 
in an average sense with in a sub cycle. 

GDPWM Algorithm: In recent years several DPWM are 
reported in [14-15] to improve the performance at high 
modulation indices. In discontinuous PWM methods 
zero sequence signals are discontinuous. During each 
sampling period, one of the phases ceases the modulation 
and the associated phase is clamped to the positive bus or 
negative bus. Hence, the switching losses of the 
associated inverter leg are eliminated. The performance 
of the PWM methods depends upon the modulation index. 
In the lower modulation range the continuous PWM 
(CPWM) methods are superior to DPWM methods, while 
in the higher modulation range the DPWM methods are 
superior to CPWM methods. However at all the operating 
modulation indices, CPWM method has higher switching 
losses than DPWM methods. Utilizing the generalized 
phase shift in the DPWM methods, generalized 
discontinuous PWM (GDPWM) is proposed in [14]. 

The SVPWM algorithm utilizing the freedom of dividing 
the zero state time within a sampling period generates 
different DPWM methods. A zero voltage vector 
distribution variable cr that divides the zero state time 
between two zero states are given by 
T 0 =T z a (5.1) 



T 7 =T z {l-a ) (5.2) 

Observation from the eight feasible switching states of 
the two level, 3-phase VSI, pole voltages V ao > Vbo > V co 
and neutral voltages V nob (k is the switching state) for 

various switching states lead to the generalized 
expressions for the pole voltages given by (6). 

1 



dc^ an 



V =— V 

v ao 2 ( 

^bo ~ V dc M bn 



Vco 2 V dc^ cn 



( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 



where M an ,M bn and M cn are the carrier based 
modulation signals comprising of fundamental frequency 
components, ranging between 1 and -1. The equations for 
the modulating signals of the top devices from (6) are 
expressed by the equation (7). 






TV- 2V 

^ r in ^ v n, 



- , i=a, b, c 



(7) 



Vdc v dc 

The neutral voltage V no can be approximated by time 
averaging over T s as given in (8). 



Vno V no \ 1 1 



T 1 +V no2 T 2 



+ v. 



no0 T 0 + V, 



no7 T 7 



= X^(T 2 - Tl ) + ^{l-2a)r z 



( 8 ) 



Expressing the switching times of the active voltage 
vectors in terms of the instantaneous maximum V max 

and minimum V min voltages, equation that is valid for all 
the six vectors, can be generalized as (9) .Using (7) and 
(9) the expressions for the modulation signals can be 
obtained [16]. The selection of cr generates infinite 
number of modulators [14], [15]. 



Vno ~b-2<j)-V min a + V max {a-l) 



(9) 



CSVPWM 



DPWM1 





DPWMMAX 



DPWM2 





DPWMMIN 



DPWM3 




a 


/ 


\ j 




/ 1 


& 


w 





Angle(radians) 



0.005 0.01 0.015 0.02 

Angle(radians) 



Figure.2.Modulation waveforms and their zero sequence signals of 
different DPWM methods 

The generalized discontinuous modulation signal is given 
by (10), 

cr = 1 - 0.5[l + sgnicos 3(cot + <?))] (10) 
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where co is the angular frequency of the reference voltage, 
8 is the modulation phase angle and is the sign 

function defined as 



sgn(x) 



+ 1, x > 0 



<0.0,x = 0 
-l,x<0 



(ID 



error in the speed is fed to the PI controller which 
generates the reference value of the torque. Based on this 
reference value, the PI controller in turn generates the 
slip speed, when added with the actual rotor speed 

generates the speed of the reference stator flux vector y/ s . 
For every sampling period y/ s and the y/* s are compared 



For CSVPWM, cr is taken as 0.5. By varying the 
modulation angle 8 , various DPWM methods can be 



generated. 



For example, 




-n 



generates 



and the reference voltage space vectors are calculated as 
given by (12). 



V ds - R,i* + 



A V d s 



( 12 . 1 ) 



DPWM1, DPWM2 and DPWM3 respectively which 

clamps for a period of 60° in every half cycle of the 
fundamental. The modulation waveforms and their zero 
sequence signals of some popular DPWM methods 
including CSVPWM methods are shown in Figure 2. 
When cr = 0 and cr = 1, 120° clamping methods 
DPWMMAX and DPWMMIN are obtained. Thus by 
varying cr and 5 the switching time periods of zero 
voltage vectors can be changed, so that different DPWM 
methods can be obtained. 

Proposed DTC: The block diagram of the proposed 
GDPWM based DTC is shown in the Figure 3. With the 
proposed method ripples in torque and flux are reduced 
significantly and with fixed switching frequency. The 
proposed DTC retains all the advantages of the basic 
DTC, such as no coordinate transformation, robust to 
machine parameter variation, no current loop etc; the 



V =R i 

v qs lv s L qs 




where, 

Ays ds -Vds ~ Vds 



AV qs V qs V qs 



( 12 . 2 ) 

(13.1) 

(13.2) 



The reference phase voltages are derived using the phase 
transformation from which imaginary switching times 
proportional to reference phase voltages are derived 
using (4). After obtaining the switching times for the 
active states and zero states the reference voltage vector 
is synthesized in an average sense with in a sub cycle. 
Thus using imaginary switching times space vector based 
GDPWM algorithm is implemented. 




6. Simulation Results and Discussion 

To verify the proposed scheme, numerical simulation has 
been carried out in Matlab/Simulink environment. 
Simulation is done using Runge-Kutta solver with a fixed 
step size of 10e-6. The average switching frequency of 
the inverter is taken as 5kHz. The reference flux is taken 
as lwb and starting torque is limited to 15 N-m. The 
induction motor used in the case study is a 1.5 KW, 1440 
rpm, 4-pole, 3 -phase induction motor having the 



following parameters: R s = 7.83Q, R r = 7.55Q, L s = 
0.475 1H, L r = 0.475 1H, L m = 0.4535 H, J = 0.06 Kg.m 2 . 
Simulation results for different DPWM based DTC 
including CDTC, showing the steady state plots are given 
in Figure. 4 - Figure. 9. Figure. 4a shows the steady state 
ripples in stator currents, speed, torque and flux with the 
CDTC. It is observed that steady state ripples in torque 
and flux are high with the CDTC. Hence to reduce the 
steady state ripples and to get the constant switching 
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frequency with reduced switching losses a GDPWM 
algorithm using imaginary switching times is proposed. 
Figures 5-10 shows the simulation results obtained from 
the proposed method. Steady state ripples and its spectral 
performance is depicted from Figureure5-Figureurel0. It 
can be observed that compared with the CDTC the 
performance in terms of ripple in torque and flux is 
significantly improved with the proposed DPWM 
methods. Also, observation of measured no load current 
waveforms and %THD’s in line current, it can be 
concluded that performance of the different DPWM 
based DTC methods is superior to the CDTC. 

Among the various DPWM methods DPWM2 is known 
for minimum switching losses. It is suitable for induction 
motor drives operating near 30° power factor angles [14]. 
Hence more elaborated simulation results showing 
different conditions of the drive, resembling transients 
during external load disturbances, starting transients, 
transients during step change in load, and transients 
during speed reversal are presented for DPWM2 based 
DTC, FigurelO. Figure 10b shows the transients in stator 
currents, speed, torque and flux during step change and 
non-linear variation in load. Figure 10c shows the 
starting transients under no-load condition. Figure lOd 
shows the transients during step change in load. Figure 
lOe shows the transients during speed reversal. Speed 
reversal command is given at 1.8sec to change the speed 
from +1300rpm to -1300rpm. 




Figure 4a. CDTC: steady state plots (Reference speed=1300 rpm). 
Fundamental (43.33Hz) = 1 .39 , THD= 21 .01 % 





Time (secs) 



Figure. 5a. DPWMMAX based DTC: steady state plots (at 1300 rpm) 
Fundamental (43.33Hz) = 1 .405 , TFID= 2.01 % 




Figure. 5b. .DPWMMAX based DTC: Measured harmonic spectra 




Time (secs) 



Figure. 6a. DPWMMIN based DTC: steady state plots (at 1300 rpm) 



Fundamental (43.33Hz) = 1 .404 , THD= 1 .99% 




Figure.6b. DPWMMIN based DTC: Measured harmonic spectra 
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Figure.7a. DPWM1 based DTC: steady state plots (at 1300 rpm) 

Fundamental (43.33Hz) = 1 .404 , THD= 1 .61 % 




Figure.7b. DPWM1 based DTC: Measured harmonic spectra 
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Figure.8a. DPWM2 based DTC: steady state plots (at 1300 rpm) 
Fundamental (43.33Hz) = 1.405 , THD= 1.65% 
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Figure.8b. DPWM2 based DTC: Measured harmonic 
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Figure 9a DPWM3 based DTC: steady state plots (at 1300 rpm) 

Fundamental (43.33Hz) = 1.404 , THD= 1.59% 




Frequency (Hz) 

Figure. 9b. DPWM3 based DTC: Measured harmonic spectra 




Figure. 10a. External load torque disturbance. 
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Figure. 10b. DPWM2 based DTC: Transients showing the external load 
torque disturbance. 
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Figure. 10d. DPWM2 based DTC: Transients during step change in load 
(A load of lON-m is applied at lsec and removed at 1.4 sec) 
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Figure. lOe. DPWM2 based DTC: Transients during reversal of speed 
(Speed reversal command is given at 1.8 sec to change the speed from 
+1300 rpm to -1300 rpm). 



7. Conclusion 

Though CDTC is simple to implement, it has few 
drawbacks like significant steady state ripples in torque 
and flux, variable switching frequency because of the use 
of hysteresis controllers. To reduce these ripples and to 
get constant switching frequency SVPWM based DTC 
came into existence. Though SVPWM based DTC could 
solve most of the tribulations, deteriorated performance 
at high modulation indices became a significant setback. 
DPWM techniques give good performance in high 
modulation regions. Also switching losses of the inverter 
are relatively not as much of Continuous PWM 
techniques. In this paper to reduce the complexity 
involved in space vector based GDPWM the concept of 
imaginary switching times is implicated. With the 
proposed PWM algorithm, infinite DPWM methods can 
be developed. Also, it can be concluded that with the 
proposed algorithm in addition to the improved 
performance at high modulation indices, switching losses 
of the inverter and computational burden involved in 
conventional GDPWM algorithm are reduced 

significantly. 
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Abstract 

We propose in this paper a Self-adaptation PI for speed 
regulation based on direct torque control (DTC) of 
Permanent Magnet Synchronous Machine (PMSM). 
Speed regulation with a conventional PI regulator 
reduces the speed control precision, increases the torque 
fluctuation, and consequentially causes low performances 
of the whole system. Using fuzzy logic method, by self- 
adaptation of conventional PI regulator parameters gives 
the appropriate Kp and Ki (proportional and integral 
coefficients respectively) which improve the system 
performance. Simulation results show that the ripples of 
both torque and flux are reduced remarkably, small 
overshooting and good dynamics of speed and torque. 

Keywords: Permanent Magnet Synchronous Machine, 
Direct Torque Control, Fuzzy PI regulator. 

1. Introduction 

The PMSM control difficulty resides in the coupling of 
control variables such as flux and electromagnetic torque. 
Two principal strategies were developed almost at the 
same time in two different research centers, Direct 
Torque Control strategy was first introduced by I. 
Takahashi, in 1986 [1]. M. Depenbrock, develop a 
similar idea in 1988 under the name of Direct Self 
Control [2]. The DTC is one of the recent researched 
control schemes based on the decoupled control of stator 
flux and torque providing a quick and robust response 
with a simple implementation in AC drives. DTC has the 
advantages of simplicity, good dynamic performance, 
and insensitive to motor parameters except the stator 
resistance. In DTC strategy, the speed sensor is not 
essential for the flux and torque estimation. Direct 
Torque Control employs two hysteresis controllers to 
regulate the stator flux and torque, which results an 
approximate decoupling between the flux and the torque 
control. The key issue of DTC scheme is how to choose a 
suitable stator voltage vector to keep the stator flux and 
torque in their hysteresis band. The conventional DTC 
disadvantages are the high torque ripples and the slow 
transient response to the step changes in torque during 



start-up. Several techniques developed to improve the 
DTC performance [3], [4], [5], [6], [7].In paper [8] the 
authors replace the conventional PI by the fuzzy logic 
regulator when the output is the torque reference. In 
paper [9] fuzzy logic is used to replace the switching 
table of DTC which make possible to choose the very 
suitable voltage vector. The paper [6] realizes Fuzzy-PI 
speed regulation for induction machine; we utilize this 
idea for PMSM in this paper. The fuzzy control is 
nonlinear and adoptive in nature, giving robust 
performances in the face of parameter variations and load 
disturbance effects. The regulator inputs are speed error 
and its change. Self-adaptation PI regulator, based on 
conventional PI regulator; consist to adjust dynamically 
the conventional PI parameters kp and ki. Simulation 
results show that the method improves the Direct Torque 
Control system performances. This paper consists of 
mathematical model of PMS machine, direct torque 
control principle; fuzzy logic technique applied to DTC, 
simulation results with interpretation, and finally 
conclusion. 



2. Motor equations in the a, p reference 
frame 

In the stationary a-ji reference frame, the model can be 
expressed as [8]: 



u a = R s i a + 



Ur =R s .ip + 



P 

Ya 

Yp 



d Ya 

dt 

d Yp 



dt 

cos 0 r - sin 6 r 
sin G r cos 6 r 



Yd 

Y q 



( 1 ) 



(2) 



Yd= L d-id + Y„ 

Yq = Lq-iq 



Ya 

Yp_ 



R d -‘a 
Lq-ip_ 



+ Ym • 



cosQ r 

sin6 r 



The mechanical equation given by: 



( 3 ) 

( 4 ) 
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^T = y(T em -T L -B m co r ) (5) 

at J 

In addition, the electromagnetic torque expressed as 
follow: 






Ya-lp-Yp-V 



0 



( 6 ) 



Where (u a , iip), (i a , ip) and (t//, r „ V.p) are the stator 
voltages, stator currents, and stator flux linkages in a, ft 
reference frame, L d , L q are the d,q axes inductance, R s is 
the stator resistance. y/ m is the flux linkage of permanent 
magnet, p is the number of pole pairs, T em is the 
electromagnetic torque, T L is the load torque, Bm is the 
damping coefficient, co r is the rotor speed and J is the 
moment of inertia. 




3. Direct torque control principles 

3.1. Flux and torque estimation and control 
In stationary reference frame, the machine stator voltage 
space vector is represented as follows: 




Where: R s , i s , y/ s stator resistance, current and flux 
respectively. V aN , V hN , V cN the three phase voltage 
inverter outputs given as follows: 

VaN=V sa =^-{2-C 1 -C 2 -C 3 ) 

V bN ^V sb = U j-(2-C 2 -C 1 -C 3 ) (9) 

V c n=V sc =^-{2-C 3 -C 2 -C 1 ) 



U c is the inverter DC supply voltage, C h C 2 , C 3 are the 
switching table outputs, and they are relevant to the 
switching strategy. 

From (7) we can estimate the stator flux as follow: 

Vs ='/,()+ I (W -tfv - 'J ( 10 ) 

So we can write: 



W sa J (y sa ^s^sa)^ 

Vsfl = J ( V s/J ~ R s i spV t 

The module of the stator flux is: 



(11) 



Ik- II = 7 



Vs 



+ Vsfl 



( 12 ) 



\f/ s is the stator flux vector, and i// s q is its initial value. 



For simplicity, it is assumed the stator voltage drop R s . is 
is small and neglected, the stator flux variation can be 
expressed as: A y/ s ~ V S .A t. 

We can estimate the electromagnetic torque using the 
following relation: 

T em =~p{Va- i p-V / p- i 3 ( 13 ) 



We can be controlled the change of torque by keeping the 
amplitude of the stator flux linkage constant and by 
controlling the rotating speed of the stator flux as fast as 
possible. We show in this section that both of the 
amplitude and rotating speed of the stator flux controlled 
by selecting the proper stator voltage vectors as shown in 
fig.l. The primary voltage vector V s , is defined by the 
equation (8) 



Figure 1. Stator flux variation in stationary (a, (3) frame 



In the DTC technique, the inverter switches obtained 
using the flux and torque errors, and the position of the 
stator flux within the six-region control of the flux. The 
flux and torque errors evaluated as follows: 

Ay/ — If/ sre f -y/ s (14) 

AT = T ref -T em (15) 



6 = tg 



-i(v± 



(16) 



Where 0 is the angle between stator flux vector and a 
axis, y/ S0[ iy/sp are lh e s t a t or flux components in (a, /3) 



reference frame. 

3.2. Switching table: 

In order to determine the inverter switching pattern using 
flux and torque errors, two hysteresis controllers are 
employing. The inverter is switched based on these errors 
and the position of the stator flux within the six-region 
control as can be seen from Table 1, in such a way that 
the inverter output voltage vector minimizes the flux and 
torque errors and determines the flux rotation direction. 



n 


1 


2 


3 


4 


5 


6 


Flux 


Couple 


K v =l 


Kcem - 1 


v 2 


V 3 


V 4 


V 5 


v 6 


Vi 


Kcem=0 


v 7 


Vo 


V 7 


Vo 


V 7 


V 0 




Kcem - “ 1 


V 6 


Vi 


V 2 
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II 
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V 7 


V 0 


V 7 


V 0 


v 7 




Kcem - “ 1 


V 5 


V 6 


Vi 


V 2 


v 3 


V 4 



Table. 1. DTC switching table. 



4. Control structure 

Figure (2) illustrates the PMSM drive scheme considered 
in this investigation. The drive consists of a Fuzzy PI 
speed controller, flux and torque controllers, space flux 
position, and PMSM. The rotor speed w r compared with 
the reference speed w re f. The resulting error and its rate of 
change are processed in the fuzzy PI speed controller for 
each sampling interval. The output of speed regulator 
considered as the reference torque T ref . 



50 



Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 8, Issue 3, ICGST, Delaware, USA, January 2009 




Figure2. A direct torque control scheme 

5. Fuzzy controller 

Speed error “Eco, ” and its rate of change AEco are using 
as inputs to the fuzzy controller. Proportional coefficient 
Kp and integral coefficient Ki are the outputs of the fuzzy 
controller. The number of fuzzy segments is chosen to 
have maximum control with a minimum number of rules. 
Triangle and trapezoidal membership functions have 
been used. The fuzzy membership functions of input and 
output variables are shown in (Fig. 3) [6]. 






Ri: If Eco is Ai, AEco is B] then KP is V t , KI is W t . 
Where A h B t , V, and ^ denote the fuzzy subsets. 
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Table. 2- Fuzzy Rules Base 



The inference method used in this paper is Mamdani’s 
procedure (inference) based on min-max decision. The 
firing strength (applied fuzzy operators) ai, for i th rules is 
given by: 

a i = mm[jU Aj ( Eco), ju B . (A Eco )) (17) 



By fuzzy reasoning, Mamdani’s minimum procedure 
gives: 



MVf ( KP ) = m in (a,- , Mv t (kp)) 
Mw t ( KI ) = min («, . MWi ( KI )) 



(18) 



Where ju A , p B , KP and KI are membership functions of 
sets A, B V and W of the variables Eco, AEco, Kp and KI, 
respectively. 

Thus, the membership functions jli v and p w of the outputs 



KP and KI are given by: 

/A-' (KP ) = niax(// v . (KP)) 

Mw, (Pf ) = max (mwMi)) 

The Maximum criterion 
defuzzification. The final 
obtained by this method. 



(19) 



method is used for 
single-valued output is 



6. Simulation results 

To verify the technique proposed, digital simulations 
based on MATLAB/SIMULINK have been implemented. 
The PMSM used for the simulations has the following 
parameters [10]: 



Figure3. The fuzzy membership functions of input and output 
variables. 



From the experience of simulation and experiment, the 
range of coefficients kp and ki are [1, 5] and [0.005, 0.02], 
respectively [6]. The speed error universe of discourse is 
divided into seven fuzzy sets. {NL (negative large), NM 
(negative middle), NS (negative small), ZE (zero), PS 
(positive small), PM (positive middle), PL (positive 
large)}, rate of change AEco includes 3 fuzzy subsets, it is 
not necessary to subdivide it, because it's changes quickly 
in DTC. Output membership KP and KI, both contain 
four fuzzy subsets as shown in (fig. 3). 

There are total of 21 rules as listed in table 2. Each 
control rule can be described using the inputs variables 
speed error “Eco”, it rate of change AEco, and output 
variables (controller parameters ki and kp). The i th rule Ri 
can be expressed as: 



Uc[V] 


350 


f[Hz] 


50 


Rs[Q] 


0.3 


Ld[mH] 


3.366 


Lq[mH] 


3.366 


J [Kg/m 2 ] 


10.8e-5 


Bm[Nm/rad/s] 


0 


\|/ m [V/rad/s] 


0.0776 


P 


5 



Table3. PMS motor parameters 



Conventional PI regulator coefficient: 
ki = 0.0553; kp =1.441. 

For proposed DTC and conventional DTC the dynamic 
responses of speed, flux, and torque for the starting 
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process without load 77=0, we applied an load torque equal to 77=7Nm at 0.6s, at t=1.7s we remove the load torque. 
We Inverse the speed at t=ls. 
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Figure4. (al) Torque response, (bl) speed response, (cl) flux 
response, for DTC with conventional PI regulator. 



Figure5. (a2) Torque response, (b2) speed response, (c2) flux 
response, for DTC with fuzzy PI regulation. 
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We compare between speed response for DTC with conventional regulator and fuzzy PI speed regulator. 



Fuzzy PI 
Conventional PI 



t(s) 



Figure6. Speed responses for DTC conventional PI regulator and DTC 
fuzzy PI regulator. 




t(s) 

Figure7. Comparison between stator flux, dashed line stator flux using 
conventional PI regulator and solid line using fuzzy PI regulator. 
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Figure8. Comparison between electromagnetic torque, dashed line 
stator flux using conventional PI regulator and solid line using fuzzy PI 
regulator. 
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Figure9. Proportional and integral coefficient estimated using fuzzy PI 
regulator. 



The simulation results show that flux and torque 
responses are very fast for two DTC methods. By 
proposed DTC technique, the ripple of torque and flux in 
steady state is reduced remarkably compared with 
conventional DTC that reduce the acoustic noise and 
vibrations. 

In fuzzy PI regulation good dynamic responses of torque 
with neglected influence of load disturbances in speed 
which restored its reference quickly. 

Figures 7 (a) and (b) represent the estimated parameters 
kp and ki of the PI regulator, we observe that kp e [1, 5] 
and ki e [0.005, 0.02], as shown in section 5. 

7. Conclusion 

In this paper, a fuzzy logic direct torque control scheme 
using fuzzy PI regulator technique is presented. Using 
fuzzy logic technique, the kp and ki can be obtained 
dynamically that gives a fast speed response. The 
simulation results suggest that FLDTC can achieve 
precise control of the stator flux and torque .Compared to 
conventional DTC, presented method the steady 
performances of ripples of both torque and flux are 
considerably improved. 
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